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ABSTRACT 

We measure the clustering of X-ray, radio, and mid-IR-selected active galactic nuclei (AGN) at 
0.2 < z < 1.2 using multi-wavelength imaging and spectroscopic redshifts from the PRIMUS and 
DEEP2 redshift surveys, covering 7 separate fields spanning ^10 square degrees. Using the cross¬ 
correlation of AGN with dense galaxy samples, we measure the clustering scale length and slope, as well 
as the bias, of AGN selected at different wavelengths. Similar to previous studies, we find that X-ray 
and radio AGN are more clustered than mid-IR-selected AGN. We further compare the clustering of 
each AGN sample with matched galaxy samples designed to have the same stellar mass, star formation 
rate, and redshift distributions as the AGN host galaxies and find no significant differences between 
their clustering properties. The observed differences in the clustering of AGN selected at different 
wavelengths can therefore be explained by the clustering differences of their host populations, which 
have different distributions in both stellar mass and star formation rate. Selection biases inherent 
in AGN selection, therefore, determine the clustering of observed AGN samples. We further find no 
significant difference between the clustering of obscured and unobscured AGN, using IRAG or WISE 
colors or X-ray hardness ratio. 


I. INTRODUGTION 

It is now well established that most galaxies host a 
supermassive black hole (SMBH) (e.g., Kormendy & 
Richstone 1995; Richstone et al. 1998; Eerrarese & Eord 
2005; Kormendy & Ho 2013). However, it is not well un¬ 
derstood what physical processes trigger intense episodes 
of accretion onto the SMBH, creating an observed active 
galactic nucleus (AGN). The broad similarities between 
the cosmic star formation history and AGN mass accre¬ 
tion history, both peaking at ^ ^ 2 and declining sharply 
at lower redshift (e.g., Soltan 1982; Madau et al. 1996; 
Eranceschini et al. 1999; Ueda et al. 2003; Zheng et al. 
2009; Serjeant et al. 2010; Aird et al. 2015), and the rela¬ 
tively tight observed correlation between SMBH mass and 
mass of the host galaxy bulge (M — a relationship; e.g., 
Magorrian et al. 1998; Gebhardt et al. 2000; Tremaine 
et al. 2002) hint at the possibility of a coeval evolution 
between SMBHs and their host galaxies. 

The vast scale difference between galaxies and SMBHs, 
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coupled with the relative rarity of the active accretion 
phase, has made it difficult to determine the physical 
mechanism(s) connecting galaxy and AGN growth. Gon- 
straining the triggering and fueling mechanism(s) of AGN 
is key to uncovering the relevant physics connecting 
SMBHs and their host galaxies. 

Clustering measurements on scales larger than a typ¬ 
ical dark matter halo (r^ > lh~^ Mpc) estimate the 
mean dark matter halo mass of AGN hosts, effectively 
placing AGN in a cosmological context (e.g.. Mo & 
White 1996; Sheth & Tormen 1999). On smaller scales 
{vp Mpc) clustering measurements estimate the 

fraction of AGN that are hosted by satellite galaxies and 
place constraints on triggering and fueling from galaxy- 
galaxy interactions and mergers. Theoretical models that 
assume different internal or external AGN triggering mech¬ 
anisms predict different large-scale clustering properties 
of AGN, as a function of both luminosity and redshift 
(e.g.. Silk & Rees 1998; Springel et al. 2005; Hopkins et al. 
2006; Hopkins & Hernquist 2009; Croton 2009; Booth 
& Schaye 2010; Eanidakis et al. 2013; Hiitsi et al. 2014). 
However, observational data suggest that only a weak lu¬ 
minosity dependence exists (e.g.. Coil et al. 2009; Krumpe 
et al. 2010; Cappelluti et al. 2010; Allevato et al. 2012; 
Koutoulidis et al. 2013). The measurement of clustering 
properties of AGN across a range of redshifts and lumi¬ 
nosities provides strong constraints to theoretical models 
of AGN. 

With the advent of the XMM-Newton and Chandra X- 
ray telescopes, early X-ray AGN clustering measurements 
at 2 ; ^ 0.5 — 2 targeted small fields and found that they re¬ 
side in massive halos from Mhaio r\j 1012-13 (Gilli 

et al. 2005; Yang et al. 2006). Later Coil et al. (2009) mea¬ 
sured the clustering of X-ray AGN sources at z ^ 1 with 
higher accuracy by using the cross-correlation of X-ray 
AGN sources with DEEP2 galaxies and using a larger field. 
They found that X-ray AGN are more strongly clustered. 
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similar to elliptical galaxies, which are more clustered 
than star-forming galaxies. Generally, X-ray AGN at 
2 ; ^ 1 — 2 are more clustered than optically-identified 
quasars as the same redshift and reside in relatively dense 
environments suggestive of being within group like envi¬ 
ronments (e.g. Gilli et al. 2005; Yang et al. 2006; Puccetti 
et al. 2006; Coil et al. 2009; Hickox et al. 2009). 

The NVSS (Condon et al. 1998) and FIRST (Becker 
et al. 1994) wide-area radio surveys identified large popula¬ 
tions of luminous, low accretion rate, mechanically-driven 
AGN (e.g., Sijacki et al. 2007). Clustering studies using 
these radio AGN found them to be strongly clustered, 
residing in very massive halos with Mhaio > 10^^ fi~^ -Ad© 
(Cress et al. 1996; Magliocchetti et al. 2004; Best et al. 
2005). Hickox et al. (2009) studied the connection between 
AGN selected using X-ray, radio, and mid-IR techniques 
by measuring the clustering, host properties, and AGN 
properties of sources in the Bootes field. They found that 
X-ray AGN and radio AGN reside in dark matter halos of 
mass Mhaio ^ Al© and Mhaio ^ Al©, 

respectively, while IR-AGN typically reside in lower mass 
halos with Mhaio < 10^^ h“^ A4©. 

The observed differences in the clustering of X-ray 
AGN, radio AGN, and IR-AGN samples indicate that it 
is crucial to test for any obscuration dependence in AGN 
clustering. The simplest unified AGN models (e.g., An- 
tonucci & Ulvestad 1985; Urry & Padovani 1995) would 
suggest that unobscured (type-1) and obscured (type-2) 
AGN should have the same distribution of environments, 
with differences in the observed obscuration due only to 
the orientation of the AGN relative to the observer. It has 
been suggested, however, that obscured and unobscured 
AGN are similar objects observed at different evolution¬ 
ary stages of SMBH accretion (e.g., Hopkins et al. 2008; 
Hickox et al. 2009). Most optical and X-ray AGN clus¬ 
tering studies do not find significant differences between 
the clustering of obscured and unobscured AGN (e.g.. 
Coil et al. 2009; Gilli et al. 2009). Hickox et al. (2011) 
found a marginal (^ 2 a) increase in the clustering ampli¬ 
tude between obscured and unobscured IR-AGN selected 
AGN at 2 ; ^ 1.25, suggesting that obscured AGN may 
reside in more massive halos. More recently, DiPompeo 
et al. (2014) and Donoso et al. (2014) found a significantly 
higher angular clustering amplitude for obscured com¬ 
pared to unobscured WISE IR-AGN at z ^ 0.9. However, 
these results measure only the angular projected cluster¬ 
ing amplitude, due to a lack of spectroscopic redshifts in 
their sample. 

Selection biases inherent in AGN identification may also 
contribute to the observed clustering signals, in that radio 
AGN are generally found in luminous, quiescent galaxies. 
X-ray AGN are found in a mixture of quiescent and star 
forming galaxies, and IR-AGN are typically found in star 
forming galaxies (e.g., Hickox et al. 2009; Aird et al. 2012; 
Mendez et al. 2013; Goulding et al. 2014). As quiescent 
galaxies are more strongly clustered than star forming 
galaxies at a given redshift (e.g., Le Fevre et al. 2005; 
Zehavi et al. 2005; Coil et al. 2008; Skibba et al. 2014), 
the observed clustering differences between AGN selected 
at different wavelengths could be due in part to differ¬ 
ences in their host populations. In order to understand 
the magnitude of this effect, one can compare the cluster¬ 
ing of AGN selected at different wavelengths to matched 
samples of inactive galaxies (e.g.. Wake et al. 2008; Man- 


delbaum et al. 2009; Coil et al. 2009; Hickox et al. 2009). 
While Coil et al. (2009) found that X-ray AGN are more 
clustered than color and magnitude matched galaxy sam¬ 
ples, Hickox et al. (2009) found that IR-AGN are less 
clustered than color and magnitude matched samples. In¬ 
terestingly, using weak leasing measurements at 2 : ^ 0.1, 
Mandelbaum et al. (2009) found that dark matter halos of 
radio-loud AGN are twice as massive as control galaxies 
of the same stellar mass and that radio AGN are more 
clustered than optically-selected AGN. 

In order to address these outstanding issues, here we 
measure the clustering properties of X-ray AGN, radio 
AGN, and IR-AGN at 0.2 < 2 : < 1.2 using the DEEP2 
and PRIMUS redshift surveys. The wealth of deep multi¬ 
wavelength data, combined with precise spectroscopic 
redshifts in these multiple fields makes this sample both 
larger and deeper than similar previous studies at these 
redshifts. We use data from multiple fields, limiting the af¬ 
fect of cosmic variance. We measure the cross-correlation 
function of AGN with dense galaxy samples, used to trace 
the large scale structure in our fields. This leads to lower 
statistical errors than measuring the auto-correlation func¬ 
tion of the AGN directly. We investigate the dependence 
of clustering with intrinsic AGN properties (e.g.. X-ray 
luminosity, specific accretion rate, hardness ratio, and 
obscuration). We create galaxy samples that are matched 
in stellar mass, star formation rate (SFR), and redshift to 
the AGN samples identified in each wavelength, to com¬ 
pare the clustering of AGN with similar inactive galaxies. 
This limits potential selection biases in comparing AGN 
samples selected at different wavelengths. 

The paper is organized as follows. In §2 we present 
the spectroscopic redshift surveys and multi-wavelength 
datasets used here. In §3 we detail the different AGN 
selection techniques and the AGN and galaxy samples 
used. In §5 we present the clustering measurements of the 
various AGN and matched galaxy samples. We discuss our 
results in §6 and conclude in §7. Throughout the paper 
we assume a standard fiat ACDM model with =0.3, 
= 0.7, and Mq = 72 km s“^ Mpc“^. 

2. DATA 

Our analysis combines multi-wavelength imaging with 
spectroscopic redshifts from the PRIMUS and DEEP2 
galaxy redshift surveys, covering eight well-known extra- 
galactic fields: the CDES-SWIRE field (Lonsdale et al. 
2003), the COSMOS field (Scoville et al. 2007), the 
DEEP2 (DEEP2; Davis et al. 2003) 02hr and 23hr fields, 
as well as the Extended Groth Strip (EGS), the Elais- 
S1 (ESI) field (Oliver et al. 2000), and the XMM-Large 
Scale Structure field (XMM-LSS; Pierre et al. 2004). We 
describe the X-ray catalogs that we use in Section 2.1, 
the radio catalogs in Section 2.2, and the mid-IR catalogs 
in Section 2.3. In Section 2.4 and Section 2.5 we briefly 
describe the PRIMUS and DEEP2 redshift surveys, re¬ 
spectively. In Section 2.6 we explain the methods used 
to estimate stellar masses and SERs for PRIMUS and 
DEEP2 sources. In Section 2.7 we provide information on 
the spatial selection function of the PRIMUS and DEEP2 
surveys that we use for our clustering analysis. 

2.1. X-ray Data 

We use existing Chandra and XMM-Newton X-ray 
source catalogs of various depths in the COSMOS, 
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Table 1 

Field information, including multi-wavelength coverage area and number of sources. 


Field 

X-ray 

Area [deg^] 
Radio IRAG 

WISE 

-^Galaxy ^ 

Number of Detected Sources 
-^Mass^ -Nx-ray^ -^Radio*^ -^Donley V\ssef 

GDFS-SWIRE 

- 

1.77 

1.77 

1.77 

20,423 

20,380 

- 

37 (41) 

131 

44 

GOSMOS 

0.93 

0.93 

0.93 

0.93 

12,284 

12,265 

203 

94 (361) 

45 

27 

Elais - South 1 

0.51 

0.90 

0.90 

0.90 

9,922 

9,903 

67 

64 433) 

59 

18 

Extended Groth Strip 

0.69 

0.71 

0.61 

0.71 

13,957 

13,178 

343 

43 (181) 

64 

15 

DEEP2 02hr 

0.58 

0.61 

0.60 

0.58 

13,222 

12,961 

61 

11 (11) 

33 

23 

DEEP2 16hr® 

0.73 

0.73 

- 

0.73 

5,645 

5,426 

31 

6 (6) 

- 

9 

DEEP2 23hr 

0.89 

0.92 

- 

0.67 

13,486 

13,239 

75 

9 (10) 

- 

19 

XMM LSS/SXDS 

2.88 

2.88 

2.84 

2.88 

35,460 

35,388 

178 

78 (141) 

157 

79 

Totals: 

7.21 

9.45 

7.64 

9.17 

124,399 

122,740 

958 

342 (894) 

489 

234 


^ Redshifts limited to 0.2 < z < 1.2 and ^iquality ^ 3. 

^ Broad-line AGN are excluded, as the optical light for these AGN is contaminated and prevents an accurate estimate of the 
stellar mass. 

X-ray detected sources with Lx > 10"^^ erg s“^. 

^ Radio AGN with Pi.4Ghz^l0^^ Watts Hz“^. Number within parentheses represents all detected radio sources with a 
robust redshift. 

® PRIMUS did not survey the DEEP2 16hr field, and we include it only for samples at 0.7 < 2; < 1.2. 


DEEP2, ESI, EGS, and XMM-LSS fields (see Aird et al. 
(2012) and Mendez et al. (2013) for details). Due to the 
large positional uncertainty of the X-ray point sources, we 
use the likelihood ratio matching technique (e.g., Suther¬ 
land & Saunders 1992; Ciliegi et al. 2003; Brusa et al. 
2007; Laird et al. 2009) to identify optical counterparts 
to each X-ray source in each field. The likelihood-ratio 
technique accounts for both the optical and X-ray po¬ 
sitional uncertainties, by calculating the probability of 
having a counterpart with a given magnitude above the 
probability of a spurious match. We place a lower limit on 
the positional uncertainty for the X-ray source location 
of 0.5" and require an optical match within 5" in any 
field. We restrict our sample to robust optical counter¬ 
parts with likelihood ratios above >0.5 and choose the 
counterpart with the largest likelihood, when there are 
multiple counterparts. Table 1 lists the area of the X-ray 
coverage in each field, as well as the number of X-ray 
sources with redshifts (see Section 2.4 for details). 

In the COSMOS field we use the public XMM-Newton 
X-ray point source catalog (Cappelluti et al. 2009; Brusa 
et al. 2010), which covers the entire 2 deg^ to a depth of 
/2-iokeV ^ 3 X 10“^^ erg s“^cm“^. We further use the 
deeper Chandra point source catalog that has a depth of 
/2-iokeV ^ 8 X 10“^^ erg s“^cm“^ and covers the central 
^0.9 deg^ (Elvis et al. 2009; Civano et al. 2012). 

In the ESI field, we use the Puccetti et al. (2006) point 
source catalog from four partially overlapping XMM- 
Newton pointings which has a depth of /2-iokeV ^ 2 x 
10-15 s“icm“^ and covers 0.52 deg^ of the PRIMUS 
area in this field. 

We use the public X-ray point source catalog from the 
deep Chandra Advanced CCD Imagining Spectrometer 
(ACIS-I) XDEEP2 survey (Coulding et al. 2012) for the 
ECS and DEEP2-02hr, DEEP2-16hr, and DEEP2-23hr 
fields. In the ECS the XDEEP2 survey contains 96 Chan¬ 
dra pointings across the field, covering an area of 0.66 
deg^. The typical full band flux limit in the merged ob¬ 
servations in this field is fx ^ 2.8 x 10“^^ erg s“^cm“^, 
though this varies across the field due to the number 
of overlapping pointings. The DEEP2-02hr, DEEP2- 
16hr, and DEEP2-23hr fields contain 12, 12 and 17 Chan¬ 
dra pointings respectively, with a full-band flux-limit of 


fx ^ 4.6 X 10 erg s ^cm ^ for all fields. In order to 
match the reported hard-band flux in the other fields, we 
convert the reported 2 — 7 keV hard X-ray band flux into 
an equivalent 2 — 10 keV hard X-ray band flux assuming 
a P = 1.9 power-law. 

In the XMM-LSS field we use the final release of the 
public XMM X-ray catalog from Chiappetti et al. (2012), 
which consists of 124 pointings of the XMM-Newton X- 
ray telescope which includes the Subaru XMM-Newton 
Deep Survey (SXDS; Ueda et al. 2008). This catalog 
contains sources to a hard-band flux limit of /x ^ 1-3 x 
10“^^ erg s“^cm“^ and fx ^ 9.3 x 10“^^ erg s“^cm“^ in 
the shallower XMM-LSS and deeper XMM-SXDS regions, 
respectively. We match the X-ray catalogs using the 
likelihood ratio matching technique described above. 

Eollowing Aird et al. (2012) and Mendez et al. (2013), 
we apply an “X-ray weight” for each X-ray detected 
source based on the ratio of the total number of X-ray 
detected sources to the predicted log(N)-log(S) relation 
of Ceorgakakis et al. (2008) at a given flux. These X-ray 
weights correct observed number densities of X-ray sources 
to the intrinsic number density and account for variations 
in the flux limit across the fields due to vignetting and 
the change in sensitivity of the telescope as a function of 
axis angle. 

2.2. Radio Data 

To select radio ACN we use existing deep Very Large 
Array (VLA) 1.4 CHz radio data in the COSMOS, EGS, 
and XMM-LSS fields. In the COSMOS field, we use the 
VLA-COSMOS Deep Project (Schinnerer et al. 2010), 
which combines the shallower data of the VLA-COSMOS 
Large Project (Schinnerer et al. 2007) with deeper cover¬ 
age in the central degree of the field. The survey provides 
radio continuum coverage for ^2,900 sources with ^ 1.5" 
resolution and a mean Icr sensitivity of 12 /iJybeam”^ 
in the central square degree and ^ 2" resolution and 
sensitivity of 15 /iJybeam”^ in the outer region. In the 
EGS, we use the AEGIS20 (Ivison et al. (2007); Willner 
et al. (2012)) VLA radio catalog which identifies 1,122 
sources from six overlapping pointings in the northern 
two-thirds of the field. The lower third of the EGS was 
not imaged due to the proximity to a bright radio source. 
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3C 295. The data were obtained from the VLA with a 5cr 
sensitivity limit of 50 /rJybeam”^ with ^ 3.8" resolution. 
In the XMM-LSS field, we use the 100 — /iJy catalog 
(Simpson et al. 2006) which contains fourteen overlapping 
pointings. The radio imaging identifies 505 radio sources 
and reaches an sensitivity limit of 12 /iJybeam”^ over 
0.8 deg^of the field. In the DEEP2-02hr, DEEP2-16hr, 
and DEEP2-23hr fields we additionally include relatively 
shallow VLA data from the Paint Images of the Radio 
Sky at Twenty-one centimeters survey (FIRST; Becker 
et al. 1995). We use the 14Mar04 catalog which contains 
946,432 radio sources above the sensitivity limit of ^ 200 
/iJybeam”^ and above the detection limit of 1 mJy. 

We use the Australian Telescope Large Area Survey 
(ATLAS) in the CDFS-SWIRE field (Norris et al. 2006) 
and ESI field (Middelberg et al. 2007). ATLAS used the 
Australian Telescope Compact Array (ACTA) at 1.4 GHz 
to survey both fields. The CDFS-SWIRE data contains 21 
pointings with 784 radio galaxies reaching a la sensitivity 
limit of ^ 40 /iJybeam“ , while the ESI data contains 
12 pointings with 1276 radio galaxies reaching a la sen¬ 
sitivity limit of ^ 30 /iJybeam”^. We find no major 
astrometric offsets between these radio catalogs and the 
PRIMUS spectroscopic catalog (described below), such 
that we assign radio counterparts to the optical redshift 
catalog by using SPHEREMATCH in IDL to identify coun¬ 
terparts within 2", corresponding to the approximately 
astrometric uncertainty in the radio catalogs. 

2.3. Mid-IR Data 

To identify mid-IR AGN, we use existing public Spitzer 
IRAC photometry in the CDFS-SWIRE, COSMOS, ECS, 
ESI, and XMM-LSS fields. IRAC provides 3.6, 4.5, 5.8, 
and 8.0 fim data which we will reference as [3.6], [4.5], 
[5.8] and [8.0]. In the CDFS-SWIRE, ESI, and XMM- 
LSS fields we use existing shallow IRAC imaging from 
Data Release 2 (DR2) from the Spitzer Wide-area In¬ 
frared Extragalactic Survey (SWIRE; Lonsdale et al. 
2003) (see Mendez et al. (2013) for details). We find no 
major astrometric offsets between these catalogs and the 
PRIMUS optical redshift catalog, and we assign IRAC 
counterparts to the optical redshift sources in all of the 
fields by matching to the closest object within 1". The 
CDFS “proper” field is not included here; instead, we 
use the larger CDFS-SWIRE field at slightly lower dec¬ 
lination, which was covered by the PRIMUS survey. In 
the COSMOS field, we reproduce the SWIRE source de¬ 
tection procedure from the SWIRE DR2 documentation 
using the IRAC mosaic images (see Mendez et al. 2013, 
for details). This ensures that we measure robust fiuxes 
and flux uncertainties using a consistent technique across 
all of our fields. For the majority of sources, our fiux 
measurements are similar to those in the S-COSMOS 
public catalog, although the public catalog tends to have 
larger uncertainties for similar brightness objects from 
the SWIRE catalogs due to their aggressive deblending 
of sources. 

In the DEEP2-02hr field we use a four-band detected 
catalog^^. The sample is drawn from Spitzer IRAC ob¬ 
servations as part of the DEEP2_CY5A/50660 program 
(PI: C. Jones). The IRAC imaging contains 34 pointings 
in each band covering the majority of the DEEP-02hr 

Catalog from A. Goulding 2013, private communication 


field. In the EGS field, we use the Barro et al. (2011) 
publicly-available IRAC [3.6] and [4.5] selected catalog. 
The catalog contains ^ 76,000 sources with [3.6] < 23.75. 
The sample is drawn from Spitzer as part of the Guar¬ 
anteed Time Observations (GTO; PI: G. Fazio) and pre¬ 
sented in Barmby et al. (2008) with additional data from 
the GO program (ID 41023; PI: K. Nandra). The GTO 
IRAC imaging comprises 52 pointing of all four IRAC 
bands over the central region of the EGS. The additional 
GO data cover the upper and lower regions of the EGS, 
banking the original strip. 

Additionally, we use data from the Wide-field Infrared 
Survey Explorer (WISE; Wright et al. 2010), which 
provides 3.4, 4.5, 12, and 22 /im photometry (bands Wl, 
W2, W3, and W4, respectively) in all of our fields. Here 
we use the public all-sky catalog from March 2012 which 
has a 5a point source sensitivity better than 0.08, 0.11, 1, 
and 6 mJy in each of the bands, respectively. We remove 
sources with spurious photometric detections and require 
sources to have SNR > 3 in the Wl and W2 bands (See 
Cutri et al. 2011, for more details). WISE surveyed 
the sky in an ecliptic polar-orbit, which increased the 
number of observations with increasing ecliptic latitude, 
causing the median coverage to vary for different fields. 
See Table 1 for the IRAG and WISE area (where at least 
Wl and W2 photometry was required) of each field. 

2.4. PRIMUS Spectroscopic Redshifts 

We use spectroscopic redshifts from the PRIMUS red¬ 
shift survey to define samples for our clustering analysis. 
PRIMUS (Coil et al. 2011; Cool et al. 2013) is the largest 
faint galaxy redshift survey completed to date, covering 
^ 9 deg^ in seven well-studied fields on the sky with multi¬ 
wavelength imaging from the X-ray to the far infrared 
(IR). The survey obtained low-resolution (A/A A ^ 40) 
spectra for ^ 300, 000 objects, targeting 80% of galaxies in 
these fields with i <22. PRIMUS used the IMACS instru¬ 
ment (Bigelow & Dressier 2003) on the Magellan-I Baade 
6.5 m telescope to observe ^ 2, 500 objects at once using 
a slitmask that covered 0.18 deg^. PRIMUS contains a 
statistically-complete sample of ^ 120, 000 spectroscopic 
redshifts to ^ab ^ 23.5. Redshifts are derived by fitting a 
large suite of galaxy, broad-line AGN, and stellar spectral 
templates to the low-resolution spectra and optical pho¬ 
tometry (see Cool et al. 2013, for details). Objects are 
classified as galaxies, broad-line AGN or stars depending 
on the best template fit. The PRIMUS redshifts are 
very precise z) ^ 0.5%) and have a low catas¬ 

trophic outlier rate, less than 3% (A^/(l z) > 0.03). 
Here we use robust (^^quaiity ^ 3, see Coil et al. (2011)) 
PRIMUS redshifts between 0.2 < z < 1.2 for the fields 
listed in Table 1. For further details of the survey design, 
targeting, and data see Coil et al. (2011); for details of 
the data reduction, redshift confidence, and completeness 
see Cool et al. (2013). 

The PRIMUS survey generally targeted all sources 
above i < 22.5 and sparse-sampled 22.5 < i < 23 sources, 
so that faint galaxy sources at the fiux limit would not 
dominate the target selection. The sampling rates are 
well defined a-priori such that building a statistically com¬ 
plete fiux-limited sample requires the tracking of both 
the “sparse sampling” weight and the “density depen¬ 
dent” weight of each object. The magnitude-dependent 
sparse sampling weight accounts for the fraction of sources 
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selected at random in the 0.5 mag interval above the tar¬ 
geting limit in each field. The density-dependent weight 
accounts for the sources in high density areas on the sky 
that are missed due to slit collisions and the finite number 
of masks observed. In these regions the observed spectra 
of adjacent galaxies would overlap on the detector if all 
galaxies were targeted (see Coil et al. 2011; Moustakas 
et al. 2013, for more details). 

Additionally, here we include a spatially-varying red- 
shift success fraction weight to account for changes in 
the observed redshift success rate across a field (i.e., due 
to differences in observing conditions for different slit- 
masks). In the PRIMUS field we use the pixelize func¬ 
tion in Mangle^^. We estimate the redshift success frac¬ 
tion by taking the ratio of highly confident sources with 
Equality > 3 to all targeted sources in the field in pixels 
of size ^ 36 arcsec^. We use a larger pixel size in the 
PRIMUS fields than in the DEEP2 fields (see Section 2.5) 
to limit Poisson noise in the shallower PRIMUS data. 

The inclusion of these targeting and completeness 
weights is important, as they correct the observations 
to more accurately represent the full galaxy population. 
However, the i ^ 23.5 selection corresponds to a rest- 
frame selection ^5000Aat ^ = 0.7, such that at the high¬ 
est redshifts we are incomplete for fainter, red galaxies. 
The redshift completeness is only a weak function of the 
color of galaxies and is a stronger function of luminosity 
(Cool et al. 2013). Eor the purposes of this paper, we 
compare the clustering of AGN samples either with each 
other or with matched galaxy samples, and these weights 
are applied to all samples. 

2.5. DEEP2 Spectroscopic Redshifts 

We also use spectroscopic redshifts from the Deep Extra- 
galactic Evolutionary Probe (DEEP2; Davis et al. 2003; 
Newman et al. 2013) redshift survey. In the DEEP2- 
02hr and DEEP2-23hr fields, PRIMUS did not target 
the 0.7 < 2: < 1.4 redshift range already covered by 
DEEP2. The combination of PRIMUS redshifts and 
DEEP2 redshifts in these fields selects galaxies uniformly 
from z = 0.2 to z = 1.4. The DEEP2 survey provides 
spectroscopic redshifts in the ECS, the DEEP2-02hr field, 
the DEEP2-16hr field, and the DEEP2-23hr field. The 
DEEP2 survey was conducted with the DEIMOS spectro¬ 
graph (Eaber et al. 2003) on the 10m Keck-II telescope. In 
the ECS, the DEEP2 survey has measured ^ 17,000 high- 
confidence redshifts (Q > 3, See Newman et al. (2013)) 
to Rab = 24.1. In the DEEP2-02hr, DEEP2-16hr and 
DEEP2-23hr fields, the survey used a photometric color se¬ 
lection to target galaxies at 0.7 < 2: < 1.4 to Rab =24.1. 
We use the Data Release 4 (DR4) catalog^^ and associ¬ 
ated window function from (Newman et al. 2013). Here 
we use redshifts between 0.2 < 2: < 1.2 in the ECS and 
redshifts between 0.7 < 2: < 1.2 in the other DEEP2 
fields. Eor all of the DEEP2 fields we require a redshift 
with a confidence greater than 95% {Q > 3). We use 
the extended optical photometry from Matthews et al. 
(2013) which contains additional Canada-Erance-Hawaii 
Telescope Legacy Survey (CEHTLS) ugriz and the Sloan 
Digital Sky Survey (SDSS) ugriz photometry matched 
to the redshift catalog. K-corrections, absolute Mb mag- 

^ ^ http://space.mit.edu/-molly/mangle/ 
http://deep.ps.uci.edu/dr4/home.html 


nitudes, and rest-frame colors are derived from K-correct 
(Blanton & Roweis 2007) from the optical photometry 
in these fields. The numbers of sources with the above 
redshift quality cuts and with estimated stellar masses 
are given in Table 1. 

We use those sources that fall within the recoverable 
spatial selection function of the DEEP2 survey. Eor the 
ECS, this precludes the use of the data from the northern 
25% of the field, which had shallower BRI photometry 
and non-uniform targeting. Eor the other DEEP2 fields 
we include all of the pointings presented in Newman et al. 
(2013). The spatial redshift success fraction reflects the 
probability that a targeted source has a secure 2:quaiity ^ 3 
redshift. Eor the DEEP2 fields we calculate this in ^ 
6 arcsec^ pixels. Using the average of six adjacent pixels 
to match the ^ 36 arcsec^ pixels used in PRIMUS does 
not change the resulting clustering measurements in these 
fields. 

2.6. iSEDfit Stellar Masses and Star Eormation Rates 

We estimate stellar masses and SERs by fitting the 
spectral energy distributions (SEDs) of our sources with 
population synthesis models using iSEDfit (Moustakas 
et al. 2013). iSEDfit is a Bayesian fitting code that 
compares the observed photometry for each source to 
a large Monte Carlo grid of SED models which span a 
wide range of stellar population parameters (e.g. age, 
metallicity, and star formation history) to estimate the 
stellar mass and SER of a galaxy. We assume a Chabrier 
(2003) initial mass function from 0.1 to 100 Al* and use 
Bruzual & Chariot (2003) stellar population synthesis 
models. We assume the following priors to construct 
the Monte Carlo grids: uniform stellar metallicity in 
the range of 0.004 < Z < 0.04; Chariot & Eall (2000) 
dust attenuation law, with an exponential distribution 
of dust, (0.25 < 7 < 2.0); an exponentially declining- 
r (0s(t) = (Al/r)e“^%) star-formation history (SEH) 
with 0.01 < r < 5.0. Stochastic bursts of star formation 
of varying amplitude, duration, and onset time are su¬ 
perimposed, allowing for a wide range of possible star 
formation histories (Kauffmann et al. 2003; Salim et al. 
2007). While a delayed-r model encompasses both a 
linearly rising {t/r <C 1) and an exponentially declining 
(t/r ^ 1) SEH history, we find no significant SER or stel¬ 
lar mass offsets or trends using different SEH models for 
our sources at 2: < 1.2, and we therefore choose to use a 
simpler model of an exponentially declining SEH. iSEDfit 
marginalizes the full posterior probability distribution of 
stellar masses and SERs over all other parameters and 
thus encapsulates both the uncertainties in the observa¬ 
tions and the model parameter degeneracies. Eor each 
source we take the median stellar mass and SER from the 
full probability distribution functions as the best estimate 
of the stellar mass and SER. The median uncertainties 
on the log stellar mass and SER are 0.08 dex and 0.2 
dex, respectively. In our analysis below we are primarily 
interested in the relative stellar mass and SER between 
sources, such that any overall offsets do not affect our 
results. 

We use iSEDfit stellar masses derived from photometry 
spanning the UV to the optical bands. Including the first 
two IRAC bands ([3.6] and [4.5]) systematically increases 
the median galaxy sample stellar mass by 0.1 dex. This is 
also the case for the X-ray detected sample; however, for 


6 


Mendez et al. 



).0 0.2 0.4 0.6 0.8 1.0 1.2 1.4 

Redshift 


26 


25 


X 

w 24 


^ 23 


22 


21 

0 , 





Radio Detected X-ray AGN 
Radio detected Donley IR-AGN 
High luminosity radio AGN 
Low luminosity radio AGN 

1.2 1.4 



IR detected Radio AGN 
Donley IR-AGN 


0.6 0.8 
Redshift 



Redshift 


Redshift 


Redshift 


Figure 1. The distributions of X-ray AGN, radio AGN, and Donley et al. IR-AGN in X-ray AGN luminosity (upper left), radio power 
(upper center), and 3.6)am IR luminosity (upper right), as well as host galaxy absolute optical magnitude {Mb) (lower left), stellar mass 
(lower center), and X-ray specific accretion rate (lower right), all as a function of redshift for 0 < 2 : < 1.4. X-ray AGN are shown with green 
circles, radio AGN with red diamonds, and Donley et al. IR-AGN with blue squares, along with inactive galaxies shown as grey dots with 
greyscale contours containing 30%, 50% and 80% of the full sample. Sources shown here have robust spectroscopic redshifts and are not 
classified as broad-line AGN. Solid vertical orange lines show the full redshift range used here (0.2 < 2 ; < 1.2), while the dashed vertical 
orange line shows the redshift used {z = 0.7) to split the full samples into higher and lower redshift samples. In the upper left panel we show 
a dashed black line for the X-ray luminosity cut used to create high and low X-ray luminosity samples, while the solid black line shows the 
lower Lx = 10^^ erg s“^ luminosity cut used for all of the X-ray AGN samples. In upper center panel the dashed black line shows the 
^i.4Ghz=10^'^ Watts Hz“^ luminosity cut that we use to create high and low luminosity radio samples. In the lower right panel we show a 
dashed black line for the specific accretion rate cut used to create high and low specific accretion rate samples. Donley et al. IR-AGN tend 
to have higher X-ray and IR luminosities compared to X-ray AGN and radio AGN sources. 


the IRAC Donley et al. IR-AGN selected sample (details 
are given below in Section 3.3) the median mass offset is 
much larger (0.5 dex). As shown in Mendez et al. (2013), 
this is due to AGN light contributing to these channels 
causing the IR-AGN to have overestimated stellar masses. 
We therefore do not include the IRAC bands in any of 
our stellar mass estimates, such that all stellar masses are 
derived using the same photometric bands, minimizing 
systematic offsets between our samples. As ^ 82% of the 
area covered by PRIMUS has GALEX UV coverage, we 
include the observed FUV and NUV photometry where 
available to improve the SFR estimates. Including the 
GALEX UV bands (compared to just using optical bands 
alone) slightly decreases the estimated stellar mass (^ 
0.02 dex) for the galaxy and AGN samples. We do 


not estimate stellar masses or SFRs for sources that are 
deemed to be broad-line AGN (BLAGN), where their 
spectra are better matched by BLAGN templates than 
by galaxy templates, as their optical photometry will be 
dominated by light from the AGN. Table 1 lists the total 
number of sources with spectroscopic redshifts in each 
field (Vgaiaxy) cind the number of sources for which we 
estimate a stellar mass (V^ass) and SFR. 

Azadi et al. (2015) investigate how the PRIMUS SFR 
estimates derived using SED fits compare with estimates 
using lOOyumdata for those galaxies detected with Her- 
schel. They find that a histogram of the SED to Herschel- 
based SFR differences peaks at zero, though there is a 
non-symmetric tail to higher SFR estimates using the 
Herschel data, such that the median difference is 0.6 dex. 
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Figure 2. Specific star-format ion rate (sSFR) versus stellar mass 
for the X-ray AGN, radio AGN, and Donley et al. IR-AGN sources. 
Golors and contours are similar to Figure 1. Dashed black lines show 
the stellar mass and sSFR cuts used to define the sSFR selected 
samples. AGN identified at different wavelengths have different 
distributions in both sSFR and stellar mass, though AGN identified 
at any wavelengths are typically identified in massive host galaxies. 
Median Icr errors on both parameters are shown in the lower left 
corner. 


Most galaxies and AGN are not Herschel-detected^ such 
that the typical difference for the full sample could be 
smaller, as Herschel detections are biased towards more 
dusty sources. However, in this paper we are concerned 
only with relative SFRs, such that any overall offsets in 
the SFR estimates will not change our results. 


we select AGN and galaxy samples using the PRIMUS 
and DEEP2 surveys in regions with either X-ray, radio or 
mid-IR imaging coverage. Below we present our selection 
criteria for our AGN and matched galaxy samples. Details 
of each AGN sample are given in Table 2. 


3.1. X-ray AGN Samples 


For the fiducial X-ray AGN sample we require that the 
detected X-ray sources have a hard-band X-ray luminosity 
Lx > 10^^ erg s~^ and a redshift in the range 0.2 < 
2 ; < 1.2. We choose to use an X-ray luminosity limit 
of Lx > 10^^ erg s“^ rather than a more conservative 
Lx > 10^^ erg s“^ limit, as this leads to larger samples 
with smaller uncertainties and no significant differences 
in our results. We have applied X-ray K-corrections ( 
(1 + F ^ 1.7 ) to estimate the hard-band X-ray 

luminosity. We create a ‘non-broadline’ subsample where 
we remove the sourced identified as BLAGN in their 
PRIMUS or DEEP2 spectra. The fiducial X-ray AGN 
sample is additionally divided into six subsamples defined 
either by an AGN property (Tx, specific accretion rate, 
or hardness ratio) or a host galaxy property (redshift, 
stellar mass, or sSFR), in order to investigate clustering 
trends with both AGN and host galaxy properties. 

For the X-ray AGN samples split by AGN luminos¬ 
ity, we divide the fiducial X-ray AGN sample into low 
luminosity {{Lx) ^ 10^^'^ erg s“^) and high luminosity 
{{Lx) ^ 10^^'^ erg s“^) samples using a luminosity cut, 
shown in the upper left panel of Figure 1. 

We also split the fiducial X-ray AGN sample by specific 
accretion rate, defined as 

AEdd = y^^ (1) 

^Edd 


_FboI_ 

1.3 X 1038 erg s-i x 0.002 ^ 


( 2 ) 


2.7. Spatial Seleetion Funetion 

In order to perform accurate clustering measurements, 
we require that all of the PRIMUS and DEEP2 sources 
used here are located within the area of each survey that 
has a well-understood spatial selection function. This 
ensures that any spatially-dependent density differences 
in the surveys that are due to target selection or missing 
data, such as in CCD chip gaps or around bright stars, 
as well accounted for In PRIMUS we require that sources 
fall within the observed window function area targeted 
with at least two slitmasks. Coil et al. (2011) provides 
details on the spatial selection function of PRIMUS, and 
Coil et al. (2004) and Newman et al. (2013) provide 
details for the DEEP2 survey. The X-ray AGN, radio 
AGN, and IR-AGN samples are identified within the areas 
with observed X-ray, radio, or mid-IR coverage. While 
there is generally overlap between the multi-wavelength 
imaging coverage, there are some areas that lack full 
multi-wavelength coverage. 

3. AGN AND GALAXY SAMPLES 

The goal of this paper is to quantify and compare 
the clustering properties of X-ray-, radio-, and mid-IR- 
selected AGN at 2 ; ^ 0.7 with each other, as well as 
with inactive galaxies with stellar mass, SFR and redshift 
distributions that match the AGN samples. To this end 


where L^dd is the Eddington limit, and LboI is the bolo- 
metric luminosity derived using the X-ray luminosity to 
bolometric luminosity relationship of Hopkins et al. (2008) 
in units of erg s“^. The specific accretion rate is a rough 
estimate of the Eddington ratio, assuming a constant scal¬ 
ing relationship between black hole mass and host stellar 
mass (e.g., AIbh ^ 0.002 Al*; Marconi & Hunt 2003). 
While there is substantial scatter in both the M — a rela¬ 
tionship and in the scaling between bulge mass and stellar 
mass of the galaxy, such that the specific accretion rate 
is not an exact estimate of Eddington ratio, it is a robust 
tracer of the rate at which the SMBH is growing relative 
to the stellar mass of the host galaxy (Aird et al. 2012). 

We create high and low specific accretion rate (A) 
samples only for AGN with a host galaxy stellar mass 
above the stellar mass limit (Al* = Mq) for which 
we are complete for quiescent galaxies at the highest 
redshifts used here and divide the fiducial X-ray sam¬ 
ple at roughly the median specific accretion rate of 
A = 10“^ erg s“^ (see the lower right panel of Eig- 

ure 1). We also create X-ray AGN samples based on 
hardness ratio, defining hard and soft samples by divid¬ 
ing the fiducial sample ad HR = 0 and requiring that 
the AGN included are identified in both the soft and 
hard X-ray bands. We adopt a simple cut of iLR = 0 as 
this approximately corresponds to Nh = 3 x 
assuming a simple absorbed power-law X-ray spectrum 










Mendez et al. 


Table 2 

Information on the X-ray, radio, and IR-AGN samples. 


AGN Selection 

Sample Name 

Redshift 

Range 

Number^ 

Density^ 

(4 

(L)- 


X-ray AGN 

Full 

0.2 - 

1.2 

958 

0.79 

0.73 

42.8 

10.7 

Radio AGN 

Full 

0.2 - 

1.2 

894 

0.55 

0.62 

23.5 

10.9 

Donley IR-AGN 

Full 

0.2 - 

1.2 

583 

0.44 

0.68 

44.2 

10.5 

X-ray AGN 

Non-broadline 

0.2 - 

1.2 

633 

0.53 

0.72 

42.7 

10.8 

Radio AGN 

Non-Broadline 

0.2 - 

1.2 

768 

0.48 

0.61 

23.5 

10.91 

Donley IR-AGN 

Non-broadline 

0.2 - 

1.2 

328 

0.25 

0.61 

44.0 

10.4 

X-ray AGN 

Low Lx 

0.2 - 

1.2 

288 

0.26 

0.57 

42.2 

10.7 

X-ray AGN 

High Lx 

0.2 - 

1.2 

328 

0.33 

0.87 

43.3 

10.8 

X-ray AGN 

Low A 

0.2 - 

1.2 

329 

0.27 

0.65 

42.4 

10.9 

X-ray AGN 

High A 

0.2 - 

1.2 

305 

0.26 

0.90 

43.4 

10.7 

X-ray AGN 

Low HR 

0.2 - 

1.2 

671 

0.54 

0.73 

42.9 

10.7 

X-ray AGN 

High HR 

0.2 - 

1.2 

287 

0.25 

0.73 

42.7 

10.8 

Radio AGN 

High Pi.4Ghz 

0.2 - 

1.2 

327 

0.21 

0.80 

24.3 

11.2 

Radio AGN 

Low Pl.4Ghz 

0.2 - 

1.2 

569 

0.34 

0.50 

23.3 

10.8 

Assef IR-AGN 

Full 

0.2 - 

1.2 

234 

0.14 

0.74 

44.5 

10.6 

Assef IR-AGN 

Obscured WISE color 

0.2 - 

1.2 

129 

0.08 

0.77 

44.5 

10.6 

Assef IR-AGN 

Unobscured WISE color 

0.2 - 

1.2 

106 

0.06 

0.70 

44.5 

10.6 


^ Number of sources in window function with applied selection cuts. 

^ Density is in units of [10“"^Mpc“^]. 

(L) is {Lx) [log(ergs“^)] for X-ray AGN samples, {PiagHz) [log(WattsH z“^)] for radio AGN, and 
{L 3 . 6 fj.Yn) [log(ergs“^)] for IR-AGN samples. 

^ Units of [log(h“^ -M©)] 


with r = 1.9at2: = 0.6 (the approximate median redshift 
of our sample), and thus roughly divides the sample into 
obscured and unobscured populations (Szokoly et al. 2004; 
Hasinger 2008). 

We further divide the fiducial X-ray AGN sample by 
various host galaxy properties, to quantify how the clus¬ 
tering of X-ray AGN depends on the host galaxy. We 
create a high and low redshift sample by dividing the 
fiducial X-ray AGN sample at ^ ^ 0.7. For both the 
stellar mass samples and sSFR samples, we require that 
the host galaxy has a stellar mass above A4* = 10^'^^ A^© 
(see Figure 2). We define the high and low stellar mass 
samples using the median stellar mass of the fiducial 
X-ray AGN sample (A4>^ ^ 10^^’^^ Al©), and we define 
high and low sSFR samples by dividing the sSFR at 
sSFR = 10.65yr“^ (see Figure 2). This cut roughly 
matches the evolving SFR-mass cut of Moustakas et al. 
(2013) at 2 ; ^ 0.7 that divides the galaxy sample into 
quiescent and star-forming galaxies. 

3.2. Radio AGN Samples 

We define four radio AGN samples based on the ob¬ 
served optical broad lines and measured radio luminosity 
of each source. For our fiducial radio AGN sample we 
require 0.2 < 2 ; < 1.2. We create a ‘non-broadline’ sub¬ 
sample where we remove the BLAGN identified by their 
PRIMUS or DEEP2 optical spectra. We have applied 
a radio K-correction ((1 -h 2 :)^^“^^; a ^ 0.5) when esti¬ 
mating the radio luminosity. Radio continuum emission 
may contain contributions from thermal bremsstrahlung 
(free-free) emission in star forming galaxies as well as 
from non-thermal synchrotron emission associated with 
radio jets emanating from an AGN. To separate these 
two populations, we follow Condon (1992) and Murphy 
et al. (2011) and define a high luminosity radio sample 
with Pi.4Ghz>10^^ Watts Hz“^, to remove any potential 
contamination from luminous starburst galaxies. Above 
this luminosity the radio emission cannot be explained 


by even extreme star formation (SER > 10^Al©yr“^) 
(Goulding et al. 2012; Hickox et al. 2009). This radio 
AGN sample reliably contains radio-loud (Class ER-II; 
Fanaroff & Riley 1974) sources. The small sample size in 
the high luminosity sample limits our analysis of the radio 
AGN sample to the fiducial redshift range (0.2 < 2 : < 1.2), 
as we do not have enough sources to create subsamples 
at different redshifts. 

This sample of radio-loud sources necessarily does not 
contain radio-quiet AGN (Mullaney et al. 2013). A va¬ 
riety of optical, mid-IR, or far-IR to radio flux ratio 
excess techniques have been suggested to identify more 
complete samples of radio-quiet AGN while limiting con¬ 
tamination from star forming galaxies (e.g. Smolcic et al. 
2008; Park et al. 2008; Donley et al. 2005; Appleton 
et al. 2004). To investigate the clustering properties of 
radio-quiet AGN. we define a low luminosity radio sam¬ 
ple (Pi.4Ghz<10^'^ Watts Hz“^). This sample includes all 
radio-detected sources below the luminosity limit, identi¬ 
fying all possible optical- or IR-excess selected sources. 

In order to investigate possible contamination of this 
low luminosity radio sample by star forming galaxies, in 
the upper panel of Figure 4 we show SFR versus radio 
luminosity for the radio detected sample. We highlight 
high luminosity (Fi.4Ghz>10^^ Watts Hz“^) sources in 
red. X-ray AGN that are radio detected in green, and 
Donley et al. IR-AGN that are radio detected in blue. 
The radio luminosity distribution of the X-ray AGN and 
Donley et al. IR-AGN samples are shown as normalized 
histograms at the bottom of the panel. Cyan points 
highlight the few highly star forming sources where their 
radio luminosity can be explained solely as due to star 
formation using the Murphy et al. (2011) SFR-to-radio 
luminosity relationship (purple line). The small number 
of sources (< 1%) above this lin shows that neither the 
high nor low luminosity radio samples are contaminated 
by star forming galaxies. The median uncertainty on the 
SFR estimates is less than 0.2 dex, suggesting that only 
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Figure 3. Stellar mass, sSFR, and redshift distributions for the X-ray AGN, radio AGN, Donley et al. IR-AGN, and galaxy tracer samples. 
The green, red, and blue filled histograms show the distributions of the X-ray AGN, radio AGN, and Donley et al. IR-AGN, respectively. 
The grey line showing the distributions for the galaxy tracer sample is scaled down by a factor of 100 for ease of comparison. Differences 
between the underlying galaxy sample and AGN identified at different wavelengths are dominated by the difficulty in selecting AGN of a 
given specific accretion rate in low mass galaxies (e.g., Aird et al. 2012). Additionally, there are substantial differences in the stellar mass, 
sSFR, and redshift distributions between AGN identified at different wavelengths. 


a few sources could be scattered above the Murphy et al. 
(2011) line. We also find that most of the X-ray and mid- 
IR AGN that are radio-detected are in the low luminosity 
radio sample Pi.4Ghz<10^^ Watts Hz“^, supporting the 
conclusion of minimal contamination from star forming 
galaxies. 

In the bottom panel of Figure 4 we show the sSFR 
versus stellar mass diagram for all PRIMUS galaxies 
(grey contours), low luminosity (black points) and high 
luminosity radio AGN (Pi .4Ghz>10^^ Watts Hz red 
points). Radio detected sources are preferentially identi¬ 
fied in massive galaxies, with the high luminosity sources 
found in the most massive galaxies. As we do not find 
large differences in the sSFRs of the low luminosity radio 
sources compared to the high luminosity radio sources, 

3.3. IR-AGN Samples 

We use the Donley et al. (2012) IRAC color-color se¬ 
lection to identify mid-IR red power law AGN. As shown 
in Mendez et al. (2013), in the PRIMUS survey this se¬ 
lection provides reliable identification of luminous AGNs 
with minimal contamination by star forming galaxies. We 
require that objects are detected in all four IRAC bands 
and have colors such that they lie within the following 
region in IRAC color-color space: 

^ = logio (, y = logic ((3) 

\j3.6fj,in/ V/4.5/Ltm/ 

x>0.08 and ^>0.15 and (4) 

^>(1.21 X x) -0.27 and (5) 

y < (1.21 X x)-b 0.27 and (6) 

/4.5^tm ^ / 3 . 6 /xm Q-Hd / 5 . 8 /xm ^ /4.5/i,m Q-Ild (7) 
fs.Ofim ^ / 5 . 8 /i m* (8) 

The small sample size of the Donley et al. IR-AGN 
sample limits our analysis to the fiducial redshift range 
(0.2 < z < 1.2). 


Additionally, we identify WISE-selected IR-AGN using 
the Assef et al. (2013) magnitude-dependent selection. 
We require sources to have measured W1 and W2 fiuxes 
such that, 

Wl-W2> 0.662 exp [0.232 {W2 - 13.97)2] (9) 

IU2< 17.11 (10) 

where W1 and W2 are in Vega magnitudes. Assef et al. 
(2013) show that this selection is 90% reliable in its identi¬ 
fication of IRAC selected AGN. This selection extends the 
Stern et al. (2012) WISE IR-AGN color selection to fainter 
limiting magnitudes, while controlling for contamination 
(see Assef et al. (2013) for details). 

In Figure 5 we compare the MIR colors of Donley et 
al. IR-AGN and Assef et al. IR-AGN. In the top panel 
we show the WISE selection plane with the Mateos et al. 
(2012) color-color selection wedge in orange and the Stern 
et al. (2012)unobscured color selection limit in red. We 
show WISE-detected Donley et al. IR-AGN with blue 
points, Assef et al. IR-AGN with red points, and PRIMUS 
galaxies with grey contours and grey outlier points. In 
this color-color space the Donley et al. IR-AGN and Assef 
et al. IR-AGN generally have similar MIR colors, with 
the Donley et al. IR-AGN extending to slightly lower 
[W1-W2] colors. In the bottom panel we show the IRAC 
color-color selection wedge with the Stern et al. (2005) 
selection in orange and the power law locus in red. We 
show Donley et al. IR-AGN with blue points, Assef et al. 
IR-AGN that are IRAC detected with red points, and the 
PRIMUS galaxy sample in grey contours. We find that 
some Assef et al. IR-AGN have high IRAC [Ch3-Ch4] 
colors, beyond the Stern et al. (2005) wedge and power 
law locus. This results from the Assef et al. IR-AGN 
selection using only the two shorter wavelength WISE 
bands and not using the longer wavelength information. 
These few sources are likely not AGN as they reside in 
the region of this diagram dominated by low-redshift 
{z < 0.3), star-forming galaxies (Mendez et al. 2013). 
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Figure 4. Host galaxy properties for radio AGN. Top: SFR 
versus radio luminosity for the fiducial radio-detected sample. High 
luminosity (-Pi.4Ghz>10^'^ Watts Hz“^) radio AGN are highlighted 
as red points. X-ray AGN and Donley et al. IR-AGN that are radio 
detected are shown with green and blue points, respectively, with 
relative distributions shown at the bottom. Radio AGN above the 
Murphy et al. (2011) SFR-to-radio luminosity relationship (purple 
line) are highlighted in cyan. In only these four sources can the 
radio light be dominated by star formation and not AGN activity. 
Bottom: Comparison of low luminosity (black) and high luminosity 
(red) radio AGN in the sSFR and stellar mass plane. The PRIMUS 
galaxy sample is shown in grey contours, with outliers beyond the 
90% contour shown with grey points. The four sources above the 
SFR-to-radio luminosity relationship of Murphy et al. (2011) are 
shown as cyan points. 

The Donley et al. IR-AGN selection uses all four bands 
to ensure a red monotonically increasing flux in the MIR, 
at the cost of requiring detections in the generally shallow 
longer wavebands. 

We divide the WISE IR-AGN sample into obscured 
and unobscured subsamples. We use the criteria of Yan 
et al. (2013), who use obscured and unobscured tem¬ 
plates at z < 1.5 to deflne a MIR-to-optical color cut of 


Figure 5. Comparison of the Donley et al. IR-AGN (blue) and 
Assef et al. IR-AGN (red) samples. Top: WISE color-color diagram 
with the Mateos et al. (2013) selection wedge shown in orange and 
the Stern et al. (2012) color selection shown in red. WISE sources 
identified by the Assef et al. IR-AGN selection technique are shown 
with red points and Donley et al. IR-AGN that are also detected in 
WISE are shown with blue points. Grey contours show all PRIMUS 
galaxies that are detected in WISE. Bottom: IRAC color-color 
diagram with the Stern et al. (2005) selection wedge shown in orange 
and the power law locus indicated in red. Assef et al. IR-AGN are 
shown with red points and Donley et al. IR-AGN are shown with 
blue points. A small number of the Assef et al. IR-AGN that are 
IRAC detected fall outside of the power law region and Stern wedge, 
due to the selection not utilizing the longer WISE wavelengths. 

(r —IE2) ^ 6 to separate these sources. Due to differences 
in the photometric Alters in our measured r-band magni¬ 
tudes, we use a synthesized SDSS r-band magnitude from 
K-Correct to ensure a uniform selection in each held. 
In the top panel of Figure 6 we show the selection of 
our samples in the optical and MIR color-color diagram, 
where obscured Assef et al. IR-AGN are shown in red and 
unobscured Assef et al. IR-AGN in blue. In the bottom 
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Figure 6. Comparison of obscured and unobscured Assef et al. 
IR-AGN. Top: Yan et al. (2013) optical and MIR color selection 
plot showing obscured (red) and unobscured (blue) Assef et al. 
IR-AGN. WISE-detected PRIMUS galaxies are shown with grey 
contours. Bottom: Optical to MIR color distribution for Assef et 
al. IR-AGN. Obscured Assef et al. IR-AGN are shown with a red 
histogram, unobscured Assef et al. IR-AGN with a blue histogram, 
and the broad line Assef et al. IR-AGN with a black histogram. 
Most broad line sources have optical-to-MIR colors of unobscured 
AGN. 

panel we show the optical to MIR color distributions of 
the obscured (red) and unobscured (blue) samples, as 
well as the optically identified broad line AGN (black) in 
these samples. We find that most broad line AGN have 
colors that identify them as unobscured AGN. 

We additionally test the Mateos et al. (2012) WISE IR 
selection technique. Similar to the Donley et al. IR-AGN 
selection, it identifies sources with a red power law in the 
mid-IR. This technique is more robust than that of Assef 
et al. (2013) as it uses longer wavelength information (W3: 
12/im) to ensure a monotonic mid-IR SED, but it is less 
complete due to the relatively shallow W3 coverage in the 


WISE survey. As we find no significant differences in the 
clustering properties of the AGN samples defined using 
Mateos et al. (2012) and Assef et al. (2013), we use the 
slightly larger Assef et al. IR-AGN sample throughout 
this paper. 

3.4. Galaxy Tracer Samples 

We use the dense galaxy samples provided by the 
PRIMUS and DEEP2 redshift surveys to measure the 
clustering of AGN using a cross-correlation measurement 
with galaxies. To do this we define galaxy “tracer” sam¬ 
ples to trace the cosmic web in the fields and at the 
redshifts of interest. Eor the fiducial galaxy tracer sample, 
we use all galaxies with robust redshifts within the fidu¬ 
cial redshift range used here, 0.2 < 2 : < 1.2. We do not 
require that the galaxy tracer sample be volume limited, 
as we are using it only to trace the large-scale structure 
in these fields; it therefore needs to span the same volume 
as our AGN samples, but it does not need to have the 
same median luminosity at all redshifts. We additionally 
split the fiducial galaxy tracer sample into low and high 
redshift subsamples for the X-ray AGN sample, split at 
the redshift of 2 : = 0.7 (Eigure 1, dashed orange line in 
each panel). This redshift cut divides the number of X-ray 
AGN into approximately equal sized samples and results 
in 30% more tracer galaxies at lower redshifts than at 
higher redshifts. 

3.5. Galaxy Matched Gontrol Samples 

We construct galaxy samples for each of the fiducial 
X-ray AGN, radio AGN, and Donley et al. IR-AGN 
samples with matched stellar mass, sSER, and redshift 
distributions. We use these as control samples to com¬ 
pare the clustering of galaxies that host AGN to similar 
galaxies without AGN. This limits the effect of AGN selec¬ 
tion biases such as AGN being preferentially identified in 
galaxies with high stellar mass (e.g., Aird et al. 2012). Ad¬ 
ditionally, comparing the clustering of AGN to matched 
galaxy samples cancels out any effective flux limits that 
arise from either the multi-wavelength AGN selection or 
the spectroscopic redshift requirement (e.g., Leauthaud 
et al. 2015). Effectively, we use individual matched galaxy 
control samples for each of the X-ray AGN, radio AGN, 
and Donley et al. IR-AGN samples to control for dif¬ 
ferences in the host galaxy properties that each AGN 
selection identifies. Coil et al. (2009) and Hickox et al. 
(2009) created similar matched galaxy samples, however 
they matched rest-frame optical magnitude and color. 
While these are easily observed properties, they are not 
as physically-motivated as stellar mass and sSER. While 
we find no significant differences in our results matching 
joint stellar mass, sSER, and redshift distributions rather 
than magnitude, color, and redshift distributions, we use 
the former parameters as they reflect intrinsic host galaxy 
physical properties. 

We measure the joint stellar mass, sSER, and redshift 
distribution for each of the X-ray AGN, radio AGN, and 
Donley et al. IR-AGN samples in stellar mass bins of 
AAI* = 0.2 dex, sSER bins of AsSFR = 0.2 dex, 
and redshift bins of Az = 0.1. The projected distribu¬ 
tions of each sample are shown in Eigure 3. Normalizing 
each distribution by the total number of galaxies in each 
bin results in the observed fraction of AGN in that bin. 
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effectively estimating the probability density in this three- 
dimensional parameter space. We use this as an estimate 
of the probability density to weight each inactive galaxy 
when creating matched galaxy samples to compare to 
various AGN samples. To ensure that we can robustly 
estimate the stellar masses and SFRs of the AGN, we 
limit the comparison of matched galaxies to AGN that 
do not have any broad-line emission in their optical spec¬ 
tra, such that the optical light is dominated by the host 
galaxy. 


4. CLUSTERING ANALYSIS 

We measure the spatial distribution of AGN using the 
two-point correlation function, which quantifies the excess 
probability above Poisson of finding two sources with a 
given physical separation. While most studies measure 
the auto-correlation function (AGE) of the AGN sample 
of interest, here we measure the cross-correlation function 
(CCE) of AGNs with galaxies, from which we then infer 
the ACF of the AGN alone. As discussed in Coil et al. 
(2009) there are two main advantages to this method. 
First, the CCF of AGN and galaxies has a much greater 
statistical power due to the larger number density of 
galaxies, which better probe the underlying dark matter 
distribution where AGN are located. Second, it does not 
require a complete understanding of the AGN selection 
function, which may not be well understood. Instead, 
all that is required is an understanding of the selection 
function of the galaxy tracer sample. 

4.1. Measuring the Cross-Correlation Funetion 

The two-point correlation function ^(r) is defined as 
the excess probability density, dP/dV^ above that of a 
Poisson random field, of a second source being physically 
separated by a distance r from a given source, 

^=n[l + C(r)] (11) 

where n is the mean number density of the sample of 
interest (Peebles 1980). The ACF measures the clustering 
of a single sample, where the two sources are from the 
same sample, while the CCF measures the clustering of 
one type of source, taken from one sample, around that 
of another type of source, taken from a second sample. 
Here we measure the CCF of AGN (A) around galaxies 
(G), which are used as a tracer sample, and find the 
excess probability above random (R) of finding AGN and 
galaxies with a given separation (r). We use the Davis & 
Peebles (1983) estimator: 




AG{r) 

AR{r) 


( 12 ) 


where AG{r) is the sum of the weighted AGN-galaxy 
pairs and AR{r) is the sum of the weighted AGN-random 
pairs, both as a function of separation. Weights are used 
to account for target selection in the PRIMUS sample 
(see Section 2.4); by applying these weights we are able to 
create a statistically-complete sample that is not subject 
to spatial biases. In the DEEP2 fields the weights are 
included in the spatial selection function which we use 
to generate the random catalogs, such that galaxies have 


unity weight. We calculate the weighted number of pairs: 




AR = Y 

^ hFAGN X A^random 


(13) 

(14) 


where icagn is the weight of a given AGN, rcgaiaxy is 
the weight of a given galaxy, Wagn is the total AGN 
weight, Wgaiaxy is the total galaxy weight, and Yrandom 
is the number of random objects. The AGN weight is 
the multiplicative combination of the targeting weight 
and any additional completeness weight such as the X-ray 
AGN weight (see Section 2.1 for details). For the DEEP2 
fields the targeting weight is unity for each source. 

Peculiar velocities distort ^(r) measurements in the 
redshift direction, along the line of sight. We therefore 
measure ^(r) in two dimensions, C(^p: 'tt), where is the 
separation perpendicular to the line of sight, which is 
unaffected by peculiar velocities, and tt is the separation 
along the line of sight. Integrating ^(r^, tt) along the tt di¬ 
mension leads to a statistic that is independent of redshift 
space distortions, the projected correlation function: 

pOO 

Wp{rp)=2 dTT^{rp,TT) (15) 

Jo 

TTmax 

dn ^{ rp , tt) (16) 

where TTmax is the maximum tt separation to which we 
integrate. As the signal to noise of C(^p: 'tt) declines 
quickly for large values of tt, we measure the projected 
correlation function by integrating to a given TTmax to limit 
shot noise. We use a larger TTmax value in the PRIMUS 
fields compared to the DEEP2 fields to account for the 
larger redshift uncertainty in the PRIMUS survey. In 
the PRIMUS fields we use TTmax = S0h~^ Mpc, while in 
DEEP2 we use TTmax = 20 Mpc. Skibba et al. (2014) 
and Coil et al. (2008) use similar values for these surveys, 
respectively. 



4.2. Jaekknife Uneertainty Estimation 

We estimate the uncertainty in our measurements us¬ 
ing jackknife resampling of the data (e.g., Lupton 1993; 
Scranton et al. 2002). For reasonably large surveys (in¬ 
cluding both PRIMUS and DEEP2) jackknife errors are 
generally similar to the cosmic variance errors in Wp de¬ 
rived from simulated mock catalogs (e.g., Zehavi et al. 
2005; Norberg et al. 2008; Coil et al. 2008; Skibba et al. 
2014). For each of our samples, we use between 10 and 
12 jackknife samples across our 8 fields, where we have 
spatially subdivided the larger fields into two or more 
subfields. The different number of jackknife samples is 
due to the multi-wave length coverage in each field (i.e. 
CDFS does not contain X-ray data; see Table I for field 
details.) We subdivide the large fields (CDFS and XMM) 
along lines of constant RA and declination such that the 
resulting subsamples probe roughly similar volumes and 
cover an area on the sky approximately equal to 1 deg^. 

The uncertainty in Wp is estimated by calculating the 
projected correlation function using each jackknife sample. 
From this collection of Wp estimates we calculate the 
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variance in the projected correlation function, 

^ '^(.u)p(rp) - Wj{rp)f, (17) 

3 

where the N is the number of jackknife samples, j indexes 
each jackknife sample, and Wj{rp) is the projected cor¬ 
relation function computed for a given jackknife sample. 
By measuring the projected correlation function using 
multiple fields across the sky, the jackknife resampling of 
fields estimates the uncertainty on our measurements due 
to cosmic variance. 


4.3. Inferring the AGN Auto-eorrelation Funetion 

Following Coil et al. (2009), we infer the AGN ACF 
from measurements of the AGN-galaxy CCF and the 
galaxy ACF. We calculate the galaxy ACF using the 
Landy & Szalay (1993) estimator: 


^{rp, tt) 


GG{rp, tt) - 2 GR{rp, n) + RR{rp, n) 

RR{rp, tt) ^ ’ 


where GG, Gi?, RR are the galaxy-galaxy, galaxy- 
random, and random-random weighted pair counts, re¬ 
spectively, where we include the galaxy targeting weights. 
We calculate the AGN-galaxy CCF and galaxy ACF in 
the same volume. We integrate the galaxy ACF projected 
correlation function to the same TTmax limits used for the 
AGN-galaxy CCF. 

We then infer the autocorrelation function of the AGN 
sample using. 


WAAiVp) 


«^ag(^p) 

WGcirp) 


(19) 


where waa is the projected AGN ACF, wag is the pro¬ 
jected AGN CCF, and wqq is the projected galaxy ACF. 
Implicit is the assumption that the spatial distributions 
of AGN and galaxies are linearly related to the underly¬ 
ing dark matter spatial distribution (i.e., that the bias 
is linear, see Section 4.5 below), and that galaxies and 
AGN are well mixed within dark matter halos. 


4.4. Power Law Fits 


The correlation function can roughly be fit by a power 
law, with f{r) = (r/ro)^, where the scale factor ro is the 
scale at which there is unity excess probability and ^ = 1. 
An analytic form can then be fit to Wp{rp): 


Wp{rp) 


rp 


/roV r(i)r(V) 

Vp) r(i) 


( 20 ) 


where F is the Gamma function. We fit this analytic func¬ 
tion to our clustering measurements in the approximately 
linear regime of = 1 — 10 Mpc. On larger scales the 
size of our fields limits the number of pair counts, which 
artificially lowers the measured correlation function and 
leads to large statistical fluctuations. 


4.5. Bias Estimation 

We use the projected correlation function to estimate 
the dark matter bias of the AGN ACF. The bias b mea¬ 
sures the relative clustering strength of the AGN sample 
to that of dark matter particles. We estimate the bias at 


the median redshift of each AGN sample using the publicly 
available code of Smith et al. (2003). We integrate the 
dark matter correlation function to a TTmax = 80 Mpc 

and then calculate the bias as 


'^^AGN 

'^^DM 


( 21 ) 


where wagn is the AGN ACF and wdm is the dark matter 
ACF on scales of = 1 — 10 Mpc. When comparing 
the clustering of different samples it is useful to compare 
the bias values instead of the clustering scale lengths, as 
the bias accounts for differences in the median redshift 
of each sample and further does not assume that ^ is a 
power law. 

Additionally, the relative bias between two AGN or 
galaxy samples is defined as the square root of the ratio 
of their respective projected correlation functions. This 
allows for a simple comparison of the clustering strength of 
two samples and is akin to comparing their absolute bias 
(relative to dark matter) values. We estimate the relative 
bias on scales of = 1 — 10 Mpc. We use the ratio of 
CCFs which does not increase the fractional uncertainty 
of the resulting bias due to the common galaxy-tracer 
ACF term in the AGN ACF. Below we present the mean 
and Icr uncertainty of the relative bias across the jackknife 
samples when comparing two samples. 


4.6. Halo Mass Estimation 

We estimate the median dark matter halo mass (Mdm) 
that hosts AGN selected at different wavelengths from the 
absolute bias measured on scales of = 1 — 10 Mpc. 
We convert the bias to the quantity u = 6c/cr{M), where 
he ^ 1.69 is the critical density for collapse and a{M) is 
the mass density fluctuation in a sphere of radius = 
(3MA/47rp) from linear theory. We use Equation (8) of 
Sheth et al. (2001) to convert the absolute bias to u and 
Equations (A8-A10) of van den Bosch (2002) to infer the 
median Mdm of the sample. 


5. RESULTS 

In this section we present the results of our correlation 
function analysis. As discussed above, we measure the 
CCF of our AGN samples relative to the galaxy tracer 
sample and the ACF of the galaxy tracer sample, in order 
to infer the ACF of the various AGN samples. We first 
present the ACF of the X-ray AGN, radio AGN, and 
Donley et al. IR-AGN in Section 5.1. We then compare 
the clustering properties of various subsamples within 
each of the X-ray AGN, radio AGN, and Donley et al. 
IR-AGN samples in Section 5.2, to investigate whether 
the clustering amplitude depends on AGN luminosity or 
hardness ratio. In Section 5.3 we present the clustering 
of obscured versus unobscured IR AGN selected in WISE. 
We compare the relative clustering between X-ray AGN, 
radio AGN, and Donley et al. IR-AGN in Section 5.4, to 
determine whether AGN selected at different wavelengths 
have different clustering amplitudes. Finally, we compare 
the clustering strength of each AGN sample relative to 
matched galaxy control samples in Section 5.6. 


5.1. Clustering of X-ray, Radio, and IR AGN 

In the left panel of Figure 7 we show the X-ray AGN 
ACF (black), which is derived from the AGN-galaxy CCF 
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Figure 7. Projected correlation functions for the X-ray AGN (left), radio AGN (center), and Donley et al. IR-AGN (right) samples. In 
each panel we show the AGN auto-correlation function (AGF; black), the galaxy-AGN cross-correlation function (GGF; red) and the galaxy 
tracer AGF (blue), as well as the uncertainties from jackknife resampling of the fields. The black line shows our fiducial sample which 
excludes the GOSMOS field, while the light grey dashed line shows the AGF including the GOSMOS field. As discussed in Section 2 the 
large over-densities in the GOSMOS field at 2 : < 1 can systematically affect the measured clustering amplitudes, as seen here. 


(red) and the galaxy tracer ACF (blue). In this figure we 
present results including (grey) and excluding (black) the 
COSMOS field. The X-ray AGN ACF is more clustered 
than the galaxy tracer sample. Table 3 lists the measured 
ro and 7 values excluding the COSMOS field, as well as 
the absolute bias and inferred median dark matter halo 
mass both with and without the COSMOS field included. 
As can be seen in Figure 7, the clustering amplitude of 
X-ray AGN increases when the COSMOS field is included; 
the bias on large scales increases by 26%. As discussed 
above, the COSMOS field contains several large structures 
at z < 1 , which both systematically increases the cluster¬ 
ing amplitude when including this field and increases the 
jackknife error. We therefore prefer to focus on results 
that exclude COSMOS when discussing the absolute bias 
of our AGN samples. The bias of X-ray AGN that we 
measure {b = 1.5 ± 0 . 2 ) corresponds to a median dark 
matter halo mass of (Mhaio) ^ h“^ A^q, which 

is generally associated with galaxy groups. The X-ray 
AGN ACF rises sharply at small projected separations 
(vp <0.7 h~^ Mpc), indicating that on small scales there 
is a increase in the number of pairs of objects within the 
same dark matter halo. 

In the center panel of Figure 7 we show the radio 
AGN ACF, again with (grey) and without (black) the 
COSMOS field included. Including COSMOS causes the 
large-scale bias to increase substantially by 50%. We find 
that similar to X-ray AGN, radio AGN are more clustered 
than the galaxy tracer sample and have a large bias value 
{b = 1.8 ± 0.1, excluding COSMOS), which suggest that 
they reside in massive dark matter halos ((Mhaio) ^ 
IO 13.3 p-i j typically associated with massive groups 
or small clusters. 

In the right panel of Figure 7 we show the Donley et al. 
IR-AGN ACF. Unlike the X-ray AGN and radio AGN, 
we find that Donley et al. IR-AGN are less clustered than 
the galaxy tracer sample. The Donley et al. IR-AGN 
sample has a lower bias value (5 = 1.3T0.5, excluding 
COSMOS) than either of the X-ray AGN or radio AGN, 
from which we estimate a median dark matter halo mass 
of ((Mhaio) 1012-8 h-^A^e). Unlike for the X-ray AGN 
and radio AGN samples, including COSMOS for Donley 
et al. IR-AGN results in a lower clustering amplitude; 


the bias decreases by 30%. 

While not shown here, we additionally measure the 
Assef et al. IR-AGN ACF and list the derived clustering 
parameters in Table 3. Similar to the Donley et al. IR- 
AGN sample, the Assef et al. IR-AGN is less clustered 
than the galaxy tracer sample and has a lower bias (6 = 
0.8 ± 0.1) than either the X-ray AGN or radio AGN. The 
large scale bias of the Assef et al. IR-AGN sample does 
not change when including the COSMOS field, though 
there is an increase in the clustering strength on small 
scales which substantially increases the slope (7 = 2 . 2 ). 
We also measured the clustering properties of Assef et 
al. IR-AGN and find similar results to both Donley et al. 
IR-AGN and Assef et al. IR-AGN. 

5.2. Relative Bias Within AGN Samples at a Given 
Wavelength 

As discussed above, the COSMOS field contains large 
over-densities at 2 ; ^ 0.3 and z ^ 0.7 (e.g., Lilly et al. 
2007; McCracken et al. 2007; Meneux et al. 2009; Kovac 
et al. 2010; de la Torre et al. 2010; Skibba et al. 2014) In 
our sample, which combines the PRIMUS and DEEP 2 
fields and therefore covers a large volume and probes a 
range of cosmic densities, we find including the COSMOS 
field systematically impacts our clustering results (some¬ 
what akin to the Sloan Great Wall, (e.g., Zehavi et al. 
2011; McBride et al. 2011)), both in terms of the mea¬ 
sured amplitude and the jackknife errors, which increase 
when COSMOS is included. 

However, including the COSMOS field is more robust 
when comparing the relative bias between two samples, as 
the same volume is used for both measurements and over¬ 
all changes in the density (i.e., cosmic variance) cancel to 
first order. This is refiected in the fact that the fractional 
jackknife errors on the relative bias values measured here 
decrease when the COSMOS field is included (which re¬ 
sults in larger samples and volumes probed). Therefore, 
when presenting relative bias measurements throughout 
the paper, we include the COSMOS field. We note that 
if the COSMOS field is excluded from our relative bias 
analysis, the significance of our results lowers (due to the 
larger errors) but the actual relative bias values do not 
change substantially. 
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Table 3 

Clustering results for X-ray, radio, and IR-AGN. 


AGN Selection 

Sample Name 



Bias‘d 


Bias^ 




without the COSMOS field 

with the COSMOS field 

X-ray AGN 

Fiducial 

4.0 ± 1.9 

1.5 ± 0.2 

1.5 ± 0.2 

12.9 

1.9 ± O.I 

13.3 

Radio AGN 

Fiducial 

5.5 ± 1.9 

1.9 ± 0.7 

1.8 ± O.I 

13.3 

2.7 ± 0.7 

13.8 

Donley IR-AGN 

Fiducial 

4.8 ± 1.3 

1.9 ± 0.3 

1.3 ± 0.5 

12.8 

0.9 ± 0.2 

II .6 

X-ray AGN 

Non-broadline 

7.0 ± 6.0 

1.4 ± 0.3 

2.5 ± 0.4 

13.6 

1.9 ± 0.2 

13.4 

Radio AGN 

Non-broadline 

5.7 ± 2.8 

1.6 ± 0.6 

1.9 ± 0.3 

13.3 

2.6 ± 0.7 

13.8 

Donley IR-AGN 

Non-broadline 

4.9 ± 1.5 

1.7 ± 0.3 

1.4 ± 0.5 

12.9 

1.1 ± 0.3 

12.5 

X-ray AGN 

Low Lx 

4.5 ± 4.3 

1.3 ± 0.3 

1.9 ± 0.3 

13.3 

2.0 ± 0.2 

13.4 

X-ray AGN 

High Lx 

4.1 ± 1.8 

1.6 ± 0.4 

1.5 ± O.I 

12.9 

1.8 ± 0.2 

13.2 

X-ray AGN 

Low A 

6.0 ± 3.6 

1.6 ± 0.4 

2.0 ± 0.2 

13.3 

2.1 ± 0.2 

13.5 

X-ray AGN 

High A 

3.8 ± I.O 

1.5 ± 0.3 

1.4 ± 0.2 

12.9 

1.8 ± 0.3 

13.3 

X-ray AGN 

Low HR 

3.2 ± 1.6 

1.4 ± 0.3 

1.4 ± 0.2 

12.8 

2.0 ± 0.0 

13.4 

X-ray AGN 

High HR 

5.0 ± 2.1 

1.7 ± 0.7 

1.6 ± 0.2 

I3.I 

1.6 ± 0.4 

I3.I 

X-ray AGN 

Low Redshift 

2.7 ± 1.8 

1.6 ± 0.4 

0.9 ± 0.2 

12.0 

1.6 ± 0.2 

13.4 

X-ray AGN 

High Redshift 

6.1 ± 2.1 

2.3 ± 0.7 

2.0 ± 0.5 

13.3 

2.7 ± 0.6 

13.6 

Radio AGN 

High Pi. 4 Ghz 

5.7 ± 2.0 

2.3 ± I.O 

1.9 ± 0.3 

13.3 

2.2 ± 0.6 

13.6 

Radio AGN 

Low Pl.4Ghz 

4.3 ± 5.0 

1.4 ± 1.2 

1.7 ± 0.3 

13.2 

2.9 ± 0.7 

13.9 

Assef IR-AGN 

Fiducial 

2.3 ± 1.3 

1.8 ± I.O 

0.8 ± O.I 

II.2 

0.8 ± O.I 

II.4 

Assef IR-AGN 

Obscured WISE color 

0.8 ± 1.6 

1.3 ± 0.5 

0.7 ± 0.2 

II.I 

0.6 ± 0.2 

II.I 

Assef IR-AGN 

Unobscured WISE color 

2.7 ± 0.9 

3.1 ± 3.8 

0.9 ± 0.2 

II.7 

1.1 ± 0.3 

12.4 


^ Correlation scale length, ro, in units of [h ^ Mpc]. 

^ Correlation power-law index 7 . 

Absolute bias estimated on scales of 1 < r^/ [h~^ Mpc] < 10 . 

Dark matter halo median mass, M^m-, in units of [log(h“^ -^o)]- 
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Figure 8. The scale dependent relative bias between various X-ray AGN samples defined by X-ray luminosity (Lx, left), specific accretion 
rate (A, center), and hardness ratio (HR, right). Error bars are derived using jackknife resampling of the fields. 




10'^ lO"" 

r [Mpc h-^ 


lO"" 10" 

r [Mpc h-M 


Figure 9. The relative bias between the high and low luminosity radio AGN samples (left), the Donley et al. IR-AGN and Assef et al. 
IR-AGN samples (center), and the obscured Assef et al. IR-AGN and unobscured Assef et al. IR-AGN (right). There is no significant 
difference between the clustering of these samples; all relative biases are consistent with unity. 


As discussed in Section 3, we divide the full X-ray AGN 
sample into subsamples depending on various AGN prop¬ 
erties: X-ray luminosity Ljy, specific accretion rate A, 
and hardness ratio. The measured clustering parameters 
for each subsample are given in Table 3. The relative bias 
between these samples is shown in Figure 8 and listed 
in Table 4 for scales associated with the one halo term 
{vp <lh~^ Mpc) and the two halo term (r^ >lh~^ Mpc). 
As above, uncertainties are derived from jackknife resam¬ 


pling of the fields, and COSMOS is included here. We 
find a 2 . 6(7 difference on small scales with Lx, such that 
the lower Lx sources are more clustered in the one-halo 
regime. For these X-ray samples defined by Lx, A, and 
hardness ratio we find no statistically significant differ¬ 
ences (> 3a) in the clustering amplitudes on either small 
or larger scales. This implies that the mass of dark matter 
host halo does not correlate with any of these properties, 
within the ranges that we probe. 
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Table 4 

Relative clustering bias for X-ray, radio, and IR-AGN samples, including the COSMOS field. 


AGN Sample Gomparison 

Relative 

0.1 < Tj 

! Bias 
. < 1 

Na 

Relative 

1 < r-p 

; Bias 
< 10 

Na 

Full X-ray AGN to Full Radio AGN ratio 

0.96 


0.04 

- 1.0 

0.87 


0.10 

-1.3 

Full X-ray AGN to Full Donley IR-AGN ratio 

1.42 


0.10 

4.4 

1.53 


0.22 

2.4 

Full Radio AGN to Full Donley IR-AGN ratio 

1.49 


0.16 

3.0 

1.81 


0.49 

1.7 

Full Assef IR-AGN to Full Donley IR-AGN ratio 

1.41 


0.24 

1.7 

1.00 


0.19 

0.0 

No Broadline X-ray AGN to Full X-ray AGN ratio 

0.95 


0.05 

- 1.0 

1.01 


0.02 

0.6 

No Broadline Radio AGN to Full Radio AGN ratio 

0.97 


0.04 

- 0.8 

0.97 


0.03 

- 1.0 

No Broadline Donley IR-AGN to Full Donley IR-AGN ratio 

1.15 


0.09 

1.6 

1.16 


0.07 

2.2 

No Broadline X-ray AGN to No Broadline Radio AGN ratio 

0.94 


0.09 

-0.7 

0.91 


0.08 

- 1.2 

No Broadline X-ray AGN to No Broadline Donley IR-AGN ratio 

1.17 


0.11 

1.6 

1.36 


0.18 

2.0 

No Broadline Radio AGN to No Broadline Donley IR-AGN ratio 

1.25 


0.20 

1.3 

1.56 


0.51 

1.1 

High Lx to Low Lx X-ray AGN ratio 

0.84 


0.06 

- 2.6 

0.95 


0.05 

- 1.0 

High A to Low A X-ray AGN ratio 

0.95 


0.18 

-0.3 

0.93 


0.09 

-0.7 

High HR to Low HR X-ray AGN ratio 

0.89 


0.15 

-0.7 

0.88 


0.09 

-1.3 

High Pi. 4 Ghz fo Low Pi. 4 Ghz Radio AGN ratio 

1.20 


0.15 

1.4 

0.87 


0.11 

- 1.2 

Obscured to Unobscured Assef IR-AGN ratio 

0.75 


0.12 

- 2.1 

0.78 


0.27 

- 0.8 

Low A X-ray AGN to No Broadline Radio AGN ratio 

1.00 


0.16 

0.0 

0.95 


0.08 

-0.7 

High A X-ray AGN to No Broadline Donley IR-AGN ratio 

1.18 


0.13 

1.4 

1.33 


0.26 

1.3 

No Broadline X-ray AGN to Matched Galaxy ratio 

0.93 


0.12 

- 0.6 

0.96 


0.05 

- 0.8 

No Broadline Radio AGN to Matched Galaxy ratio 

0.91 


0.17 

-0.5 

1.09 


0.07 

1.3 

No Broadline Donley IR-AGN to Matched Galaxy ratio 

0.89 


0.08 

-1.5 

0.80 


0.16 

-1.3 


In the left panel of Figure 9 we show the relative bias 
between high and low luminosity radio AGN. We do not 
find significant differences in the clustering of high and 
low luminosity radio AGN on either small or large scales 
(see Table 4). The lack of a dependence of the clustering 
amplitude on radio luminosity may be surprising, given 
that the highest luminosity radio AGN are found in the 
most massive quiescent galaxies compared to lower lumi¬ 
nosity radio AGN (bottom panel of Figure 4). While the 
lower luminosity radio AGN sample contains more galax¬ 
ies with slightly lower mass (A4* ^ AI©, compared 
to ^ Mq for the high luminosity sample), the 
clustering signal is dominated by the most massive ob¬ 
jects in the sample. Additionally, we find no significant 
difference in the sSFR distribution of the host galaxies 
as a function of the radio luminosity, which suggests that 
high and low luminosity radio AGN have fairly similar 
host galaxies. 

For each of the AGN properties tested above (i.e., lu¬ 
minosity, specific accretion rate) we test the significance 
of the clustering differences when comparing sources in 
the upper and lower quartiles, as opposed to the upper 
and lower halves. This allows us to look for clustering 
differences between the most extreme sources in each 
parameter of interest. This does not change any of our 
results, however, the larger uncertainties that result from 
using smaller samples may prevent us from detecting any 
underlying differences. 

In the center panel of Figure 9 we show the relative bias 
between the Donley et al. IR-AGN and Assef et al. IR- 
AGN samples. We list the measured clustering parameters 
in Table 3 and relative bias between these samples in 
Table 4. We find that these samples have consistent 
clustering properties, given the error bars, however the 
bias of Donley et al. IR-AGN is 37% higher than that 



Redshift 

Figure 10. Redshift distributions for the obscured (red hatched), 
unobscured (blue), and broad line (dark grey hatched) Assef et al. 
IR-AGN. The obscured and unobscured IR-AGN populations have 
substantially different redshift distributions; the obscured sample 
peaks at 2 : ~ 0.7, while the unobscured sample has a much flatter 
distribution, including more sources are low redshift. The higher 
redshift sources in both samples are more likely to have broad 
optical lines. 


of the Assef et al. IR-AGN. Both samples have lower 
estimates for the median dark matter halo mass than the 
X-ray AGN or radio AGN. 

5.3. Relative Bias of Ohseured versus Unohseured 
WISE IR-AGN 
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Following Donoso et al. (2014) and DiPompeo et al. 
(2014), we compare the clustering of obscured and un¬ 
obscured WISE-selected AGN. In the right panel of Fig¬ 
ure 9 we show the scale dependent relative bias between 
obscured and unobscured Assef et al. IR-AGN. The 
obscured Assef et al. IR-AGN are less clustered than 
the unobscured Assef et al. IR-AGN at the ^ 2cr level 
on small scales and at the ^ la level on larger scales. 
Within the uncertainties, therefore, we do not find a sig¬ 
nificant difference in the clustering amplitudes of these 
samples. This suggests that the differences found by 
Donoso et al. (2014), where the COSMOS field alone was 
used to determine the redshift distributions of these two 
populations, was impacted by using a single field. Here, 
using eight fields that cover ^ 9.1 deg^ of the sky, and 
using spectroscopic redshifts for each source, we do not 
find a significant difference in their clustering. 

To understand this further, we show in Figure 10 the 
redshift distributions of the obscured (red), unobscured 
(blue), and broad line Assef et al. IR-AGN (black) in 
our sample. While the median redshift of the obscured 
Assef et al. IR-AGN ((z) ^ 0.77) and unobscured Assef 
et al. IR-AGN ((z) ^ 0.70) are similar, the samples 
have very different redshift distributions. The obscured 
AGN peak at (z) ^ 0.7, whereas the unobscured AGN 
have a flatter distribution and peak at both lower and 
higher redshift, at (z) ^ 0.3 and (z) ^ 1.0. Within the 
unobscured sample, the broad line sources typically have 
higher redshifts than the non-broad line sources. 

The differences in the redshift distributions of obscured 
and unobscured Assef et al. IR-AGN strongly limits 
the interpretation of angular clustering measurements of 
these two populations. While there are broad similarities 
between our redshift distribution and that of Donoso 
et al. (2014), our obscured sample shows a much broader 
redshift distribution extending out to (z) ^ 1.2. Our 
unobscured sample also has a larger fraction of sources 
at low redshift ((z) ^ 0.3). Our redshift distribution is 
more similar to the redshift distributions found in the 
Bootes Survey (DiPompeo et al. 2014). 

Similarly, we test both the Donley et al. IR-AGN sam¬ 
ple and Assef et al. IR-AGN samples for any dependence 
of the clustering amplitude with obscuration and find no 
significant differences. For Donley et al. IR-AGN we test 
both an optical-to-WISE color cut {r — W2 = 6.0) and 
an optical-to-IRAC color cut (r — [4.6] =6.1) and find no 
significant differences using either cut. 

5.4. Comparison of the Clustering of X-ray, Radio, and 

IR AGN 

In Figure 11 we show the relative bias between AGN 
identified at different wavelengths; the results are given in 
Table 4. We find that radio AGN are more clustered than 
X-ray AGN (red line) on large scales (15% higher bias), 
but the difference is not significant (1.3cr). Comparing the 
X-ray AGN and Donley et al. IR-AGN samples (purple 
line), there is a significant difference (4.4 a) on small 
scales, where X-ray AGN have a 38% higher bias, while 
on large scales X-ray AGN having a 44% higher bias 
(2.4cr). Comparing the radio AGN and Donley et al. IR- 
AGN samples (cyan line), radio AGN have a 40% higher 
bias on small scales (3.0cr) and an 58% higher bias on 
large scales (1.7cr). The relative bias averaged over all 
scales results in a difference at the 2.5cr level. 



Figure 11. The relative bias between the X-ray, radio, and Donley 
et al. IR-AGN samples with error bars determined from jackknife 
resampling of the fields. We show the relative bias between X-ray 
to radio AGN (dark red). X-ray to Donley et al. IR-AGN (purple), 
and radio to Donley et al. IR-AGN (cyan). We find that the X-ray 
and radio AGN have similar clustering amplitudes, both of which 
are larger than that of Donley et al. IR-AGN. 

On small scales (r^ <lh~^ Mpc) the correlation func¬ 
tions of both the X-ray AGN and radio AGN are signif¬ 
icantly higher than that of the Donley et al. IR-AGN. 
This may be due in part to the difficulty in identifying 
Donley et al. IR-AGN in quiescent galaxies, due to the 
1.4 jam stellar bump entering into the mid-IR photometry 
(Mendez et al. 2013). This selection effect limits the num¬ 
ber of Donley et al. IR-AGN that can be identified in 
quiescent host galaxies, which would decrease the cluster¬ 
ing amplitude on all scales, though particularly on small 
scales (reflecting differences in color-dependent clustering; 
e.g., Zehavi et al. 2005; Coil et al. 2009; Skibba et al. 
2014). 

5.5. Dependence on Specific Accretion Rate 

As shown in Hickox et al. (2009) and Mendez et al. 
(2013), there is a substantial difference in the specific ac¬ 
cretion rate (A) distributions of AGN selected at different 
wavelengths. To account for these differences in A, here 
we compare radio AGN and Donley et al. IR-AGN to 
X-ray AGN with similar A values. In the large left panel 
of Figure 12 we show A versus stellar mass for X-ray AGN 
(green circles). X-ray AGN that are also radio AGN are 
shown with red diamonds, and those that are also Donley 
et al. IR-AGN are shown with blue squares. We divide 
X-ray AGN into a high A sample and a low A sample at 
A = 10“^ erg s“^ AIq . Most Donley et al. IR-AGN lie 
above this line, while most radio AGN lie below this line. 
Here we do not include sources below At* ^ At©, 
to ensure that we are roughly complete in stellar mass at 
all redshifts. We also remove all broad-line sources from 
this comparison, as we require a stellar mass estimate for 
the high-A X-ray AGN and low-A X-ray AGN samples. 

In the upper-right panel of Figure 12 we show the scale 
dependent relative bias between the high-A X-ray AGN 
and Donley et al. IR-AGN. In the lower-right panel of 
Figure 12 we show the relative bias between the low-A 
X-ray AGN (green) and radio AGN (red). We find no 
significant differences between either sample. The limited 
significance of these results are dominated by the low 
number of non-broadline AGN in each sample. 

5.6. Comparison with Matched Galaxy Control Samples 

As discussed above, AGN samples identified at different 
wavelengths are biased in terms of identifying specific 
types of AGN in specific types of host galaxies. In general. 
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Figure 12. Dependence of clustering amplitude on specific accretion rate (A) for X-ray, radio, and Donley et al. IR-AGN. Left panel: 
Specific accretion rate versus stellar mass for X-ray detected AGN. X-ray AGN that are also radio AGN are shown with red diamonds, and 
those that are also Donley et al. IR-AGN are shown with blue squares. The X-ray luminosity limit for our sample is shown as a light-grey 
dashed diagonal line. X-ray AGN below the Al* = 10^'^^ Al© stellar mass limit (vertical dashed line) are shown with open symbols, while 
X-ray AGN with broad lines in their optical spectra are shown as small grey points. The A = —2 line defines the high A and low A X-ray 
AGN. Top right panel: Relative bias between the high A X-ray AGN and Donley et al. IR-AGN for the non-broadline samples. Bottom 
right panel: Relative bias between the low A X-ray AGN and radio AGN for the non-broadline samples. We find that Donley et al. 
IR-AGN are similarly clustered as high A X-ray AGN while radio AGN are clustered similarly as low A X-ray AGN. 


AGN are more easily identified in more massive galaxies 
(e.g., Silverman et al. 2011; Aird et al. 2012). Additionally, 
there can be substantial differences in the sSFR and 
redshift distributions of the host galaxies of AGN selected 
at X-ray, radio, and IR wavelengths (e.g., Hickox et al. 
2009; Coil et al. 2009; Mendez et al. 2013; Goulding et al. 
2014). 

Differences in the host galaxy populations can influence 
the observed AGN clustering amplitude, which must be 
understood before AGN clustering can be used to con¬ 
strain the AGN triggering mechanism. To account for 
this, we compare the clustering of each of our AGN sam¬ 
ples to that of matched galaxy control samples that have 
the same stellar mass, sSFR, and redshift distributions 
as the AGN samples. 

As this requires robust estimates of the stellar masses 
and SFRs, we remove optical broad-line AGN for these 
comparisons, in order to limit AGN contamination of 
the optical broadband photometry used in the SED fits. 
While broad line AGN are a substantial fraction (34%) 
of the X-ray AGN population, excluding them does not 
significantly change the measured clustering properties. 
As seen in Table 3, excluding broad line AGN from the 
fiducial X-ray AGN sample leads to a 1% change in the 
large scale bias. This implies that at least for the Lx 
range and redshift range considered here, the clustering of 
narrow line and broad line X-ray AGN are not significantly 
different. The radio AGN sample contains the smallest 
fraction of broad line sources (14%). Excluding the broad 
line radio sources results in only a 3% change in the bias 
on large scales. The broad line fraction of the Donley et al. 
IR-AGN sample is also substantial (31%), similar to that 
of the X-ray AGN sample, and excluding these sources 
leads to a 14% change in the bias. The relative biases 
between the X-ray, radio, and IR AGN samples after the 
broad line sources have been removed do not change too 
substantially, but the fractional errors do increase, due 
to the smaller sample sizes. 



Figure 13. The relative bias between the X-ray AGN (green), 
radio AGN (red), and Donley et al. IR-AGN (blue) and their 
respective stellar mass, sSFR, and redshift matched galaxy control 
samples. We do not find any significant differences in the clustering 
properties of any of the AGN and matched galaxy samples; all 
relative biases are consistent with unity. 

In Figure 13 we show the relative biases between (non- 
broadline) AGN samples identified at different wave¬ 
lengths and their matched control galaxy samples; the 
results are listed in Table 4. We find no significant dif¬ 
ferences in the clustering amplitude of either the X-ray 
AGN, radio AGN, or Donley et al. IR-AGN and their 
matched galaxy control samples on small or large scales; 
all differences are significant at 1.5cr or less. 

This suggests that the physical effect (s) of the host 
galaxy large-scale environment is either sub-dominant in 
AGN triggering or is not separable from the host galaxy 
properties. The strong agreement between the cluster¬ 
ing of AGN host galaxies and similarly-selected inactive 
galaxies, on both small and large scales corresponding 
to the one- and two-halo terms, indicates that the AGN 
triggering mechanism either does not act on these scales 
or correlates with the properties of the galaxies in which 
they are identified. It could be possible, for example, that 
an environmental effect that triggers AGN also causes 
changes in the host galaxy’s sSFR, such that active and 
inactive galaxies with the same sSFR distribution have 
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the same clustering properties. However, as shown else¬ 
where, there are not substantial correlations between host 
galaxy sSFR and AGN incidence, once selection effects 
are taken into account (e.g., Aird et al. 2012 ; Azadi et al. 
2015). 

6 . DISCUSSION 

We have combined spectroscopic redshifts with multi¬ 
wavelength imaging to quantify the clustering properties 
of AGN selected at different wavelengths and compared 
their clustering with matched galaxy samples. In this 
section, we discuss the implications of these findings. 
In Section 6.1 we compare our results to other multi¬ 
wavelength AGN clustering studies in the literature. In 
Section 5.5 we investigate whether AGN clustering ampli¬ 
tude depends on specific accretion rate, and in Section 6.2 
we discuss the halo mass dependence of AGN activity. 

6 . 1 . Comparison with the literature 

In Figure 14 we compare the measured clustering am¬ 
plitude of AGN identified at different wavelengths using 
the bias parameter estimated on scales lh~^ Mpc < rp < 
10 h~^ Mpc. We do not compare ro and 7 values in the lit¬ 
erature, due to the degeneracy between these parameters. 
In the top left panel we compare the bias parameters for 
our three full samples selected at different wavelengths, 
while each of the other panels compares our results at a 
given wavelength with other results from the literature. 
In all panels the filled, colored symbols show our results 
excluding the COSMOS field, while the open, colored 
symbols include the COSMOS field. 

In the top right panel we compare the bias values of var¬ 
ious X-ray AGN clustering studies. For each of the results 
shown that include the COSMOS field (our open symbol 
samples, Gilli et al. (2009), Allevato et al. ( 2011 ), and 
Koutoulidis et al. (2013)), the median dark matter halo 
mass of X-ray AGN is above (Mhaio) ^ 10^^h“^ Al©, 
while those that do not include the COSMOS field (Gilli 
et al. (2005), Yang et al. (2006), Krumpe et al. (2010), 
Hickox et al. (2009), and Coil et al. (2009)) find a median 
mass below (Mhaio) ^ 10^^-^ h“^ Ad©. Both Gilli et al. 
(2005) and Yang et al. (2006) use one or two small fields 
(< 0.5 deg^) to measure the clustering of X-ray AGN, 
which leads to an underestimate of the cosmic variance 
in their measurements. Including the COSMOS field sys¬ 
tematically raises the bias due to the large over-densities 
found in the field at 2 : < 1 . Gilli et al. (2009) find an 
increase of ^ 24% in the bias when including the over¬ 
density found at 2 : ^ 0.36 within the COSMOS field. We 
find a similar ^ 20 % increase. When we exclude COS¬ 
MOS, we find similar results to Coil et al. (2009) and 
Hickox et al. (2009). 

Interestingly, Allevato et al. (2011) and Koutoulidis 
et al. (2013) include the COSMOS field in their results 
and find X-ray AGN bias values similar to our fiducial 
results, which exclude COSMOS. However, Allevato et al. 
( 2011 ) derive their results using only the COSMOS field 
and bootstrap errors within that field, such that their 
errors due to cosmic variance are underestimated. The 
analysis in Koutoulidis et al. (2013) spans a very wide 
redshift range, 0 < 2: < 3, and the results shown here are 
for the median redshift of their sample. This makes a 
direct comparison with results derived in smaller redshift 
bins somewhat difficult. 


Taken together, all of these results show that X-ray 
AGN are typically found in somewhat more massive dark 
matter halos at {z) ^ 0.9 compared to {z) ^ 0.4. The 
lower redshift X-ray AGN have a lower median X-ray lu¬ 
minosity {{Lx) ^ 10 ^^’^ erg s“^) compared to the higher 
redshift X-ray AGN {{Lx) ^ 10^^’^ erg s“^). However, 
we do not find a correlation between clustering ampli¬ 
tude and X-ray luminosity in our samples, which sug¬ 
gests that the luminosity differences between the red¬ 
shift samples is not driving the difference in clustering 
strength. We also do not find a significant difference in 
the median stellar mass of the lower redshift AGN hosts 
in our sample ((Ad*) ^ 10 ^^’^ Ad©) compared to the 
higher redshift hosts ((Ad,,) ^ 10^°-^ Ad©). This differ¬ 
ence is similar to that of the low X-ray luminosity AGN 
((Ad*) ^ Ad©) and high X-ray luminosity AGN 

((Ad*) ^ -^o) samples, suggesting that differences 

in stellar mass are not driving the redshift-dependent re¬ 
sults seen here. 

In the bottom left panel we show the bias of radio 
AGN compared to results in the literature. Generally, the 
bias we measure agrees well with other published studies 
and indicates a relatively high dark matter halo mass of 
(Mhaio) ^ h“^ Ad©. Including the COSMOS field 

again increases the bias by ^ 50% but also substantially 
increases the error bars, due to that one field having a 
systematically different clustering amplitude compared to 
the other fields. The Wake et al. (2008) results are higher 
than other measurements and are derived from the cross¬ 
correlation of radio-loud luminous red galaxies (LRGs) 
(A. 4 Ghz> 10 ^^ Watts Hz“^) with the main 2 SLAQ LRG 
survey. Their radio-detected LRGs have luminosities of 
3-51/*, far higher than the average radio AGN, which may 
account for their high bias value. Our low luminosity 
radio AGN have a consistent bias value as the Donoso 
et al. (2010) and Hickox et al. (2009) samples. 

In the bottom right panel we show the IR-AGN selected 
using the Donley et al. IR-AGN selection technique (blue 
diamond) and the Assef et al. IR-AGN selection technique 
(red square). We additionally show the Assef et al. IR- 
AGN obscured (red diamond) and unobscured (red circle) 
subsamples. We compare our results to those of Hickox 
et al. (2009, grey diamond, ^ ^ 0.5) and Gilli et al. 
(2007, grey square, ^ > 0.7), as well as angular clustering 
estimates from Donoso et al. (2014) and DiPompeo et al. 
(2014) for obscured and unobscured sources. While the 
individual selection techniques compared (e.g. Donley 
et al. IR-AGN, Stern et al. IR-AGN, / 24 /xm-selected) 
differ, we generally find that the IR-AGN samples have 
lower bias values than the X-ray AGN or radio AGN and 
therefore typically reside in lower mass dark matter halos 
((Mhaio) ^ 10^^-^ h“^ Ad©). We do not find a significant 
difference between the clustering amplitudes of samples 
identified using the Donley et al. IR-AGN or Assef et al. 
IR-AGN techniques. This is in agreement with Hickox 
et al. ( 2011 ), who consider somewhat higher redshift 
((2:) ^ 1 . 2 ) and higher luminosity {{Lx) ^ 10^^ erg s“^) 
IR-AGN than we do here. However, our results coupled 
with theirs suggest that there is not a significant difference 
in the clustering of obscured and unobscured IR-AGN for 
a range of different redshifts and luminosities. 

We also find no significant difference in the bias of the 
obscured and unobscured Assef et al. IR-AGN. While 
there is a ^ 29% higher bias for the unobscured sources 
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Figure 14. Comparison of the absolute bias as a function of redshift for AGNs identified at different wavelengths in the literature. Top 
left panel: Bias of the X-ray AGN (green circles), radio AGN (red squares), and Donley et al. IR-AGN (blue diamonds) presented in this 
paper, where the uncertainties are derived using jackknife resampling of the fields. In all panels, we show the estimated bias for samples 
in this paper. The filled symbols show the measured values excluding the COSMOS field. Black lines show constant dark matter halo 
mass (Mhaio = 10^^ At© - M^aio = h“^ -M©). Top right panel: Comparison of bias values in the literature for X-ray AGN. 

Our results are shown in green, while grey points show results from the literature to 2 ; ~ 1. Bottom left panel: Bias of the high luminosity 
(red squares) and low luminosity radio AGN samples (red circles). Grey points show results from the literature for radio AGN. Bottom 
right panel: Bias of the Donley et al. IR-AGN (blue diamond) and Assef et al. IR-AGN samples (red square). We show both the Assef et 
al. IR-AGN WISE obscured subsample (red thin diamond) and unobscured subsample (red thin circle). Grey points show results from the 
literature for IR-AGN. 


compared to the obscured sources, this is not significant 
Icr). Our results do not agree with the angular cluster¬ 
ing measurements from DiPompeo et al. (2014) or Donoso 
et ah (2014). Since we use spectroscopic redshifts, our 
results are more robust to differences in the redshift dis¬ 
tributions of the two samples, which could be driving the 
difference in the clustering amplitude of the obscured and 
unobscured samples in these other papers, as discussed 
above. Additionally, we find that the unobscured Assef et 


al. IR-AGN have a brighter median W1 flux (^ 0.5dex) 
than obscured Assef et al. IR-AGN, which suggest that 
these samples have different effective luminosities, which 
will result in different redshift distributions. Since we 
require spectroscopic redshifts our samples are smaller 
and therefore our statistical error bars are larger, however 
the systematic errors associated with our spectroscopic 
samples should be much lower. 

We note that the redshift success rate (the fraction of 
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Figure 15. The relative bias for galaxy samples matched in stellar 
mass, sSFR, and redshift to AGN samples selected at different 
wavelengths. Error bars are from jackknife resampling of all fields. 
We show the scale-dependent relative bias for the X-ray AGN to 
radio AGN matched galaxies (dark red), X-ray AGN to Donley et 
al. IR-AGN matched galaxies (purple), and radio AGN to Donley 
et al. IR-AGN matched galaxies (cyan). 

PRIMUS targets for which we derive a robust redshift) 
is very similar for the WISE IR-AGN sample 72%) as 
for the full PRIMUS sample 75% Cool et al. (2013)), 
when we account for the number of WISE IR-AGN sources 
outside the redshift range of PRIMUS, using the DiPom- 
peo et al. (2014) WISE IR-AGN redshift distribution 
from the AGES survey (the AGN and Galaxy Evolution 
Survey Kochanek et al. 2012). We find a small trend 
(10% difference from the median) between the redshift 
success fraction and the observed W2 magnitude, where 
brighter WISE IR-AGN have a higher redshift success 
fraction. We find a similar trend for both the obscured 
and unobscured sources, which would bias both samples 
to somewhat more luminous sources. 

6.2. Does AGN activity depend on halo mass? 

We show that differences in the host galaxy popula¬ 
tions of AGN identified at different wavelengths likely 
contribute to (if not fully account for) the differences 
in the inferred host dark matter halo masses of these 
AGN. To account for the known host galaxy selection 
biases in AGN identification, we compare the clustering 
of each AGN sample to galaxy control samples with the 
same distribution of stellar mass, sSER, and redshift in 
Section 5.6 and find no significant differences. Likewise, 
Leauthaud et al. (2015) find that X-ray AGN do not 
reside in different halos than inactive galaxies, when they 
control for the stellar mass distribution of their detected 
X-ray AGN. The lack of significant differences in the 
clustering amplitudes of AGN and their galaxy control 
samples at even the < 2 a level strongly suggests that the 
physical mechanisms that are fueling and triggering AGN 
either correlate with galaxy environment on scales much 
smaller than those that we probe here (r^ ^ 0.1 Mpc) 
or that AGN triggering is not correlated with the mass 
of the dark matter halo. Eor example, our results are 
consistent with Ellison et al. (2011), who use pairs of 
optically-selected AGN to identify a sharp increase in 
the activity of AGN at close separations (< 80 h~^ kpc); 
here we are sensitive to > 100 h~^ kpc only. 

The validity of these statements relies on our ability to 
estimate the relevant galaxy properties and the fractional 
uncertainty in the clustering measurements. Both the 
stellar mass and sSER estimates that we use improve 
upon previous techniques found in the literature, as we fit 
full stellar population synthesis models to the broad-band 


SEDs of the AGN and galaxy samples in a consistent 
manner, to limit biases that result from using individual 
bands or simple color-to-SER correlations (e.g., Mostek 
et al. 2013). Additionally, we use relatively large samples 
of AGN and galaxy control samples with spectroscopic 
redshifts and maximize their measurement power by cross- 
correlating each with the PRIMUS and DEEP2 galaxy 
samples. This lowers the fractional uncertainty in our 
measurements and increases the significance of our results. 

While we do not find a significant difference between 
the clustering of radio AGN and matched galaxy control 
samples. Wake et al. (2008) and Donoso et al. (2010) find 
that radio AGN are more clustered than their stellar mass- 
matched samples. Both compare mass-matched luminous 
red galaxies (LRGs) to radio AGN at 2 : ^ 0.5 and find 
that radio AGN are significantly more clustered than the 
matched sample. While both authors examine radio AGN 
with stellar masses and halo masses higher than those 
probed by our sample, neither of these papers explicitly 
control for differences in the sSER distributions, which we 
find to be important in comparing different AGN samples 
to their host galaxies. 

Einally, the matched galaxy control samples account 
for the individual selection biases from the AGN samples 
identified at different wavelengths. These biases depend 
on the depth of the sample and the different wavelengths 
that are used to identify the AGN, both of which lead to 
differences in host galaxy properties (e.g., stellar mass, 
sSER, and redshift). In Eigure 13 we show the scale 
dependent relative bias between matched galaxy control 
samples for our AGN identified at different wavelengths. 
Similar to the relative biases found between the AGN 
samples, we find that the clustering strength of both the X- 
ray AGN and radio AGN matched galaxy control samples 
are higher than the Donley et al. IR-AGN matched galaxy 
control sample on all scales. The clustering strength of 
the radio AGN matched galaxy control sample is higher 
than the X-ray AGN matched galaxy control sample on 
small scales but is statistically similar on large scales. 

The consistency between the matched galaxy control 
samples and the AGN identified at different wavelengths 
can be used to better measure the clustering of AGN. Dif¬ 
ferences in the host galaxy properties of AGN selected at 
different wavelengths can be understood entirely as being 
due to selection effects (i.e., Mendez et al. 2013). We have 
shown here that these selection biases can entirely ac¬ 
count for differences in the observed clustering properties 
of AGN selected at different wavelengths. This confirms 
and extends the X-ray AGN results of Leauthaud et al. 
(2015) to AGN detected in the radio and IR as well. The 
clustering of AGN can therefore be understood in terms 
of the clustering of their host galaxy populations. 

7. GONGLUSIONS 

In this paper we measure the clustering properties of 
X-ray, radio, and IR-selected AGN in the PRIMUS and 
DEEP2 spectroscopic surveys. Within the X-ray AGN 
sample we measure the dependence of clustering on X- 
ray luminosity, specific accretion rate, and hardness ratio. 
Within the radio AGN sample we measure the dependence 
of clustering on radio luminosity, and within the IR- 
AGN sample we measure the dependence of clustering on 
obscuration. As the AGN in these samples span a wide 
range of specific accretion rates (which roughly traces 
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Eddington ratio), we also investigate the dependence of 
clustering on specific accretion rate. We quantify the 
relative clustering strength (or relative bias) between 
each AGN sample, as well as between the AGN samples 
and galaxy control samples that are matched in stellar 
mass, sSFR, and redshift. The main results from our 
work are as follows: 

1. The clustering amplitude of observed X-ray AGN, 
radio AGN and Donley et al. IR-AGN at 0.2 < ^ < 
1.2 differ, indicating that they reside in different 
mass dark matter halos. X-ray AGN and radio 
AGN cluster similarly, and both are more clustered 
(at ^2cr) than Donley et al. IR-AGN, especially 
on scales < 1.0 Mpc. We estimate that our 
X-ray AGN, radio AGN, and Donley et al. IR- 
AGN samples have median dark matter halo masses 

of Mhalo ~ 1012-9 \C\Mq, Mhalo - 1013-3 h”! Mq, 

and Mhalo ^ 10^^’^ h Ad©, respectively. 

2. We find no significant dependencies (< 2cr) on the 
clustering amplitude with X-ray luminosity, spe¬ 
cific accretion rate, or hardness ratio. We also 
find no significant difference in the clustering ampli¬ 
tude of radio-loud AGN (Pi.4Ghz>10^^ Watts Hz“^) 
compared to low luminosity radio-detected AGN 
(^1.4Ghz<10^^ Watts Hz“^). 

3. There is no significant difference in the clustering of 
IR-AGN samples selected using either the Donley 
et al. (2012) or Assef et al. (2013) selection tech¬ 
niques. Using either selection we find no significant 
difference in the clustering amplitude of obscured 
versus unobscured IR-AGN, using WISE-optical 
colors to define obscuration. This suggests that 
previously published differences determined using 
angular clustering are dominated by differences and 
uncertainties in the redshift distributions of these 
sources. 

4. The clustering amplitudes of the X-ray AGN, radio 
AGN and Donley et al. IR-AGN samples agree well 
with those of the matched galaxy control samples, 
which have the same distribution in stellar mass, 
sSFR, and redshift of the AGN host galaxy samples. 

It is now understood that all AGN selection techniques 
have inherent biases no matter which waveband or tech¬ 
nique is used. For example, AGN identified using X-ray, 
radio, or MIR emission, as used here, are all more eas¬ 
ily detected in massive host galaxies. This means that 
clustering results should always be interpreted as the 
clustering of observed AGN samples, above a given flux 
limit and therefore corresponding to a given stellar mass 
limit. It is clear that in addition to selection biases with 
stellar mass, there are additional biases with respect to 
the SFR of the host galaxy, where radio AGN tend to 
be identified in quiescent galaxies and IR-AGN have a 
bias towards being detected in star-forming host galaxies. 
This affects the observed clustering of the AGN, which 
should only be interpreted relative to matched galaxy 
samples. The full population of AGN amongst galaxies 
of all stellar masses is likely to exhibit different clustering 
properties that the observed AGN, which are more easily 
identified in massive galaxies. When we match samples 


with respect to stellar mass, SFR, and redshift we find 
excellent agreement between the clustering of AGN host 
galaxies and inactive galaxies. Therefore AGN clustering 
can be understood entirely in terms of galaxy clustering 
(and the dependence of clustering on galaxy properties) 
and AGN selection effects. Looking forward, theoreti¬ 
cal models that include AGN evolution and predict the 
clustering of AGN must include AGN selection biases in 
order to accurately constrain the physical mechanisms 
triggering AGN. Future observational results from surveys 
such as the Dark Energy Survey (DES) and the Hyper 
Suprime-Cam (HSC) Survey will also need to account for 
these AGN selection effects as they push to lower AGN 
luminosities in order to tighten constraints on theoretical 
models. 
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